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^ ! Abstract 

'NT 

•/^ i We consider the effect of the coupled variations of fundamental constants on the nucleon mag- 

netic moment. The nucleon g-factor enters into the interpretation of the measurements of variations 
in the fine-structure constant, a, in both the laboratory (through atomic clock measurements) and in 

r^ I astrophysical systems (e.g. through measurements of the 21 cm transitions). A null result can be 

translated into a limit on the variation of a set of fundamental constants, that is usually reduced to a. 
However, in specific models, particularly unification models, changes in a are always accompanied 
by corresponding changes in other fundamental quantities such as the QCD scale, Aqcd- This work 
tracks the changes in the nucleon g-factors induced from changes in Aqcd and the light quark masses. 

^^ . In principle, these coupled variations can improve the bounds on the variation of a by an order of 

5— ( I magnitude from existing atomic clock and astrophysical measurements. Unfortunately, the calcula- 

tion of the dependence of g-factors on fundamental parameters is notoriously model-dependent. 

1 Introduction 

Any definitive measurement of a temporal or spatial variation in a fundamental constant, such as the 
fine-structure constant a, would signal physics beyond the standard model, and in particular a violation 
of the equivalence principle which is one of the foundations of general relativity. In many cases, such an 
observation would indicate the existence of a new light (usually scalar) degree of freedom [1]. Indeed, 
there has been considerable excitement during the last decade over the possible time variations in a from 
observations of quasar absorption systems [2, 3, 4, 5, 6, 7]. 

In effectively all unification models of non-gravitational interactions, and certainly in models in 
which one imposes gauge coupling unification at some high energy scale, a variation in a is invari- 
ably accompanied by variations in other gauge couplings [8, 9]. In particular, variations in the strong 
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gauge coupling, a^, will induce variations in the QCD scale, Aqcd, as can be seen from the low energy 
expression for Aqcd when mass thresholds are included 



, mc mb mt \ 27 
Aqcd = /" ( 3 ) exp 
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27r 
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(1) 



for a renormalization scale ^ > ?7it up to the unification scale [8, 9, 10], where r7T,c,b,t are the masses of 
the charm, bottom, and top quarks. Because fermion masses are proportional to hv where ft, is a Yukawa 
coupling and v is the Higgs vacuum expectation value (vev), variations in Yukawa couplings will also 
affect variations in Aqcd so that 

AAqcd „Aa 2 f Av A/ic A/ib AhA 

*^'^° - R h— 3 — + — -^ + — -^ + . (2) 
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Typical values for R are of order 30 in many grand unified theories, but there is considerable model- 
dependence in this coefficient [11]. 

Furthermore, in theories in which the electroweak scale is derived by dimensional transmutation, 
changes in the Yukawa couplings (particularly the top Yukawa) lead to exponentially large changes in 
the Higgs vev. In such theories, the Higgs expectation value is related to the Planck mass, Mp, by [9] 

/ 27rc\ 
w~Mpexp( j, (3) 

where c is a constant of order 1, and at = /it /4vr. For c ~ /it ~ 1, 

~^^— , (4) 

V /it 

with S ~ 160, though there is considerable model-dependence in this value as well. For example, in 
supersymmetric models, S can be related to the sensitivity of the Z gauge boson mass to the top Yukawa, 
and may take values anywhere from about 80 to 500 [12]. This dependence gets translated into a variation 
in all low energy particle masses [13]. 

In addition, in many string theories, all gauge and Yukawa couplings are determined by the expecta- 
tion value of a dilaton and we might expect [9] 

Ah _lAa 
h la 

assuming that all Yukawa couplings vary similarly, so that they all reduce to h. Therefore, once we allow 
a to vary, virtually all masses and couplings are expected to vary as well, typically much more strongly 
than the variation induced by the Coulomb interaction alone. 

Irrespective of the purported observations of a time variation in a, many experiments and analyses 
have led to limits on possible variations [14, 15]. Furthermore, the use of coupled variations has led to 
significantly improved constraints in a wide range of environments ranging from big bang nucleosyn- 
thesis [9, 16, 17, 18, 19, 20, 21], the Oklo reactor [22, 24], meteoritic data [22, 23, 24], the microwave 
background [25, 20] and stellar evolution [26]. 

This article explores the possibility that the strongest existing limits on the fine-structure constant, 
namely those derived from atomic clock measurements, can also be enhanced by considering such cou- 
pled variations. We expect the effect of induced variations in Aqcd and the light quark masses to enter 
through the nucleon magnetic moment. Existing experimental limits on a from atomic clock experi- 
ments assume constant ;Up „. Indeed, limits on the variations of quark masses in units of the QCD scale, 
i.e. rriq/AQCD, from atomic clock measurements have been derived [27, 28]. Given a (model-dependent) 
calculation of the nucleon magnetic moment (or equivalently its (7-f actor), we can derive sharper bounds 
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on the variation of a from existing data. Unfortunately, because of the model-dependence, we find that 
while the limits are generally improved (by as much as an order of magnitude), there is considerable 
uncertainty in the precise numerical limit. As a corollary, we apply our results to astrophysical measure- 
ments such as those which rely on the 21 cm line which also depends on /ip,n- 

The article is organized as follows: In section 2, we outline the procedure of obtaining limits on 
a from atomic clock experiments. In particular, we examine the detailed dependence on the nuclear 
^-factors which will be subject to variation. In section 3, we derive the dependence of the nucleon mag- 
netic moment on Aqcd and the light quark masses. Because there is no unique (or rigorous) method for 
calculating baryon magnetic moments, we consider several different approaches. The most straightfor- 
ward employs the constituent quark model. Surprisingly, this model is quite effective in matching the 
observed baryon magnetic moments. Even within this broad approach, our result will depend on the 
calculation of the nucleon mass, as well as the calculation of the constituent quark mass; each carrying a 
significant degree of uncertainty. We also consider an approach based on chiral perturbation theory, and 
a method based partially on lattice results. In section 4, we apply these results to atomic clock measure- 
ments and derive "improved" limits on the variation of a. Finally, in section 5, we extend these results 
to measurements involving the 21 cm line and summarize our results. 

2 Atomic clock constraints 

2.1 From frequency shifts to constants 

The comparison of atomic clocks provides a constraint on the relative shift of the frequencies of the two 
clocks as a function of time, on time scales of the order of a couple of years. This observation (or lack 
thereof) can be translated into a constraint on the time variation of a fundamental constant. Using QED, 
the frequency of the atomic transitions can be expressed (see e.g. [29]) in terms of the fine structure 
constant a, the electron-to-proton mass ratio, /i = rric/mp and the gyromagnetic factor gi = 2^j/|Un» 
where Hi is the nuclear magnetic moment, and ^un = 2^ i^ '^^e nuclear magneton. 
The hyperfine frequency in a given electronic state of an alkali-like atom is given by 

J^his - RooC X Ams X giXa^ X fix F^sia), (6) 

where Rq^ the Rydberg constant, Ahfs is a numerical factor depending on the atomic species and Fhfs(a) 
is a factor taking into account relativistic corrections (including the Casimir contribution) which depends 
on the atom. We omitted the effect of the finite nuclear radius on hyperfine frequency in Eq. (6), since 
the effect of varying the nuclear radius is shown to be smaller [30, 31] than the effects of varying other 
parameters which we consider in this work. Similarly, the frequency of an electronic transition is well- 
approximated by 

t'elec ^ RooC X ylelec X FekciZ, a), (7) 

where, as above, A^iec is a numerical factor depending on each particular atom and Feiec is the function 
accounting for relativistic effects, spin-orbit couplings and many -body effects. Even though an electronic 
transition should also include a contribution from the hyperfine interaction, it is generally only a small 
fraction of the transition energy and thus should not carry any significant sensitivity to a variation of the 
fundamental constants. 

Relativistic corrections are important [32] and are computed by means of relativistic A^-body cal- 
culations [33, 34, 35, 36]. These can be characterized by introducing the sensitivity of the relativistic 
factors to a variation of a defined by 

SlnF 

Ka = -J-^ • (8) 

ma 
The values of these coefficients for the transitions that we shall consider below are summarized in Table I . 
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Table 1: Sensitivity of various transitions on a variation of the fine structure constant. From Refs. [33, 



34, 35, 36]. 



Atom 


Transition 




Sensitivity Kq 


^H 


Is -2s 




0.00 


87Rb 


hf 




0.34 


133Cs 


25l/2(^ = 2) - (F . 


= 3) 


0.83 


171 Yb + 


^^1/2 - ^^3/2 




0.9 


199Hg + 


'S,/2 2 AV2 




-3.2 


87Sr 


5*0 — ' -fo 




0.06 


27^1 + 


5*0 — ' -Fb 




0.008 



2.2 Experimental constraints 

Over the past several years, many comparisons of atomic clocks have been performed. We consider only 
the latest result of each type of comparison for our analysis. 

• Rubidium: The comparison of the hyperfine frequencies of rubidium and caesium in their elec- 
tronic ground state between 1998 and 2004 [29] yields 



±^j!^\ 



. In ( ^-^ ) = (0.5 ± 5.3) X 10^16 yr~^ 



(9) 



From Eq. (6), and using the values of the sensitivities Ka, we deduce that this comparison con 
strains 

i^Cs ^Cs 0.49 



oc a 

i^Rb 9Rh 



(10) 



Atomic hydrogen: The Is — 2s transition in atomic hydrogen was compared to the ground state 
hyperfine splitting of caesium [37] in 1999 and 2003, setting an upper limit on the variation of vu 
of (—29 lb 57) Hz within 44 months. This can be translated in a relative drift 



± In (^] = (32 ± 63) X W~'' yT-\ 
dt V ^H / 



(11) 



Since the relativistic correction for the atomic hydrogen transition nearly vanishes, we have z^h 
i?oo so that 

(12) 



^Cs 2.83 

cc gcsjJ-a 



Mercury: The i99jjg+ '^Si/2 — ^D^/2 optical transition has a high sensitivity to a (see Table 1) so 
that it is well suited to test its variation. The frequency of the ^^^Hg^ electric quadrupole transition 
at 282 nm was thus compared to the ground state hyperfine transition of caesium first during a two 
year period [38] and then over a 6 year period [39] to get 






-3.7 ±3.9) X lO^^^yr^^ 



(13) 



While z^cs is still given by Eq. (6), j^jg is given by Eq. (7). Using the sensitivities of Table 1, we 
conclude that this comparison test the stability of 



'^Cs 6.03 
0C5rcs^a • 



(14) 
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• Ytterbium: The "^81/2 — ^^3/2 electric quadrupole transition at 688 THz of ^^^Yb+ was compared 
to the ground state hyperfine transition of caesium. The constraint of [40] was updated, after a 
comparison over a six year period, which leads to [41] 

— Inf^^l = (0.78 ±1.40) X 10"^^ yr-^ (15) 

This tests the stability of 

^Cs 1.93 .,^. 

• Strontium: The comparison of the ^^o — ^Pq transition in neutral ^^Sr with a caesium clock was 
performed in three independent laboratories. The combination of these three experiments [42] 
leads to the constraint 

Alnf ^ ) = (1.0±1.8) X lO^i^yr-^ (17) 

dt V '^Sr / 

Similarly, this tests the stability of 

^Cs 2.77 /ION 

ocfifcs/ua . (18) 

• Atomic dyprosium: The electric dipole (El) transition between two nearly degenerate opposite- 
parity states in atomic dyprosium should be highly sensitive to the variation of a [34, 35, 43, 44]. 
The frequencies of two isotopes of dyprosium were monitored over a 8 months period [45] showing 
that the frequency variation of the 3.1-MHz transition in ^^^Dy and the 235-MHz transition in 
^^^Dy are 9.0ib6.7 Hz/yr and -0.6ib6.5 Hz/yr, respectively. This provides the constraint 

" (-2.7±2.6) X 10"^^yr-\ (19) 



a 



at lo" level, without any assumptions on the constancy of other fundamental constants. 

Aluminium and mercury single-ion optical clocks: The comparison of the ^Sq — ^Pq transition in 
^^Al+ and "^8112 — ^^5/2 i^i ^^^Hg+ over a year allowed one to set the constraint [46] 

d i„ /^ '^Al A f r o , ^7 r^^ , , -, r^-17 -1 



_ In ^ii = (-5.3 ± 7.9) X 10~" yr'^ (20) 



Proceeding as previously, this tests the stability of 



fHg 

which, using Eq. (21) directly sets the constraint 

a 



""^^ oc a3.208^ (21) 



-1.65 ± 2.46) X lO"^'^ yr~\ (22) 

a 



since it depends only on a. 



Experiments with diatomic molecules, as first pointed out by Thomson [47] provide a test of the 
variation of /i. The energy difference between two adjacent rotational levels in a diatomic molecule is 
inversely proportional to Mr^"^, r being the bond length and M the reduced mass, and the vibrational 
transition of the same molecule has, in first approximation, a \/M dependence. For molecular hydrogen 
M = ?Ttp/2 so that the comparison of an observed vibro-rotational spectrum with a laboratory spectrum 
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gives an information on the variation of rup and 771^. Comparing pure rotational transitions with elec- 
tronic transitions gives a measurement of /U. It follows that the frequency of vibro-rotation transitions is, 
in the Born-Oppenheimer approximation, of the form 

I'^Ej (c^;^^ + c^.^ ^ + c^„, n) , (23) 

where c , , c .. and c , are some numerical coefficients. 

dec ' vib rot 

The comparison of the vibro-rotational transition in the molecule SF6 was compared to a caesium 
clock over a two-year period, leading to the constraint [48] 

A In (^^] = (-1.9 ± 0.12 ± 2.7) x 10"^^ yr-\ (24) 

where the second error takes into account uncontrolled systematics. Now, using Table 1 again and Eq. (6) 
for Cs, we deduce that for a vibrational transition, 



l^SF6 



^gcs^/J2a^-^\ (25) 



2.3 Nuclear ^r-factors 



All the constraints involve only 4 quantities, fi, a and the two gyromagnetic factors gcs and g^\^. It 
follows that we need to relate the nuclear ^-factors that appeared in the constraints of the previous sub- 
section, with the proton and neutron gi-factors that will be calculated in Section 3. 

An approximate calculation of the nuclear magnetic moment is possible in the shell model and is 
relatively simple for even-odd (or odd-even) nuclei where the nuclear magnetic moment is determined 
by the unpaired nucleon. For a single nucleon, in a particular {l,j) state within the nucleus, we can write 

/ ^^31 +9s , J j = / + 1 

«={5il[2(/ + lte-..l "" {,=1-1 <''> 

where gi = 1(0) and gs = gp{gn) for a valence proton (neutron). 

From the previous discussion, the only (/-factors that are needed are those for ^^Rb and ^^^Cs. For 
both isotopes, we have an unpaired valence proton. For ^^Rb, the ground state is in a P3/2 state so that 
/ = 1 and i = |, while for ^'^^Cs, the ground state is in a gj/2 state corresponding to / = 4 and i = |. 
Using Eq. (26), the nuclear (/-factor can easily be expressed in terms oi gp alone. Using gp = 5.586, we 
find g = 7.586 for ^^Rb and g = 3.433 for ^^^Cs, while the experimental values are g = 5.502 for 
^^Rb and g = 5.164 for ^^^Cs. 

The differences between the shell model predicted (/-factors and the experimental values can be 
attributed to the effects of the polarization of the non- valence nucleons and spin-spin interaction [27, 31]. 
Taking these effects into account, the refined formula relevant for our discussion of ^^Rb and ^^^Cs is 

g = 2 [(/n b {s,y + ((/p - 1)(1 - 6)(s.)° + j] , (27) 

where (/„ = —3.826, (s^)" is the spin expectation value of the single valence proton in the shell model 
and it is one half of the coefficient of gs in Eq. (26), and b is determined by the spin-spin interaction and 
it appears in the expressions for the spin expectation value of the valence proton (sz^) = (1 — b){sz)° 
and non-valence neutrons (sz^) = b{sz)°- Following the preferred method in [27, 31], it is found 

and 

{s-.,) = {szy-{s-z:). (29) 
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Therefore, the variation of the g'-factor can be written as 



9 gp g gn g b g 



From Eq. (28), (29) and (30), we find, by using the experimental (/-factors, 

5p ' 9n 



0.764^^ - 0.172^^^ - 0.379-^ 



5Rb 
Sgcs _ 
gcs 

2.4 Summary of the constraints 



.6b 
~b 



-0.619^^ + 0.152-^ + 0.335— 





9p 



gn 



(30) 

(31) 
(32) 



Given the discussion in the two previous subsections, and in particular Eqs. (31) and (32), the atomic 
clock experiments give constraints on the set {gp,gn, b, /i, a} and thus variations in the relative frequency 
shift UAB = i'a/^b are given by 



6VAB 
l^AB 



A 



9p 



5p 



5gn 

9n 



5b 
J 



+ \— + A^ h Aq 



6n 



6a 
a 



or equivalently 



l^AB 



Ao„ h Ao„ h Afe- + A^— + Aq — , 



'9p" 



with the coefficients {A^p, A^,^, Af,, A^, A^} summarized in Table 2. 
For the sake of comparison, the shell model gives 

^ ^ 0.736^ 



a 



(33) 



(34) 



and 



5Rb 



9p 



-1.266 



5ffp 



5Cs 9'p 

The main difference arises from the dependence in g^ and b but the order of magnitude is similar. 



(35) 



(36) 



Table 2: Summary of the constraints of the atomic clock experiments and values of the coefficients 
{\gp,\g^,Xb, A^, Xa} entering the decomposition (34). 



J'ab/j^ab (y^ 



Clocks 



J^AB 



A, 



■9jiL. 



A 



yn 



Xb 



X, 



A. 



(0.5 ±5.3) X 10"^*' 

(32 ± 63) X 10-16 
(-3.7 ±3.9) X 10-16 
(0.78 ± 1.40) X 10-1^ 

(1.0 ±1.8) X 10-15 

(-1.9 ±0.12 ±2.7) X 10-1^ 

(-2.7 ±2.6) X 10-15 

(5.3 ± 7.9) X 10-1'^ 



9Cs ^,U.49 



2.83 



Cs-Rb 
Cs-H 

CS - 199Hg+ 

Cs - i'^iYb+ 
Cs-Sr 
Cs - SFg 
Dy 

199Hg+ . 27A1H 



9R.b 

gcs ^ a 
gcs ^ a' 

5Cs ^ a' 



,6.03 
1.93 

2.77 
2.83 



a 



a 

-3.208 



-1.383 
-0.619 
-0.619 
-0.619 
-0.619 
-0.619 







0.325 
0.152 
0.152 
0.152 
0.152 
0.152 







0.714 
0.335 
0.335 
0.335 
0.335 
0.335 








1 
1 
1 
1 

0.5 





0.49 
2.83 
6.03 
1.93 

2.77 
2.83 

1 
-3.20J 
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3 Nucleon magnetic moments, current quark masses and Aqcd 

In this section, we will review several approaches in the literature in calculating the nucleon magnetic mo- 
ments, including the non-relativistic constituent quark model (NQM), chiral perturbation theory (xPT), 
and a method combining the results of ^PT and lattice QCD. We will try to extract the dependence of the 
nucleon magnetic moments on the current quark masses and Aqcd from the expressions given by each 
of these approaches. 

3.1 The non-relativistic constituent quark model approach 

The NQM, which approximates hadrons as bound states of their constituent quarks gives a good approx- 
imation to the measured baryon magnetic moments [49]. In this model, the baryon magnetic moments 
are expressed in terms of the Dirac magnetic moments of their constituent quarks, with the coefficients 
given by the baryon spin/flavor wave functions. For the proton and neutron, the magnetic moments are 

4 1 4 1 

fJ'P = o^^u - o/"d and /in = g/Wd - g/Wu , (37) 

where Hu = ^JIT ^^'^ ^^<i ~ "32^- Here, M^ and Mj are the constituent u and d quark masses, 
respectively, with their values around a third of the nucleon mass, to be compared with the much smaller 
u and d current quark masses, ttIu and m^, which are several MeV. For the three light flavors (n, d and 
s), the main part of their constituent quark masses have a strong interaction origin, with the dynamics of 
the virtual gluons and quark- antiquark sea being responsible for the large masses [50], while the current 
quark masses which contribute only a small portion of their corresponding constituent quark masses are 
of pure electroweak origin. 

From Eq. (37), the nucleon magnetic moment in units of the nuclear magneton /xn = 2m~^ ^^^ ^^' 
the gi-factor of the nucleon, can be written as 

(TTl TTl \ 

where c^ = 8/9 and cj = 1/9 for the proton, and c^ = —2/9 and Cd = —4/9 for the neutron. In 
the study of hadron properties, the constituent quark masses are usually taken as fitting parameters, 
with Mu = Md often assumed [49], since isospin is a good approximate symmetry. We will assume 
this relation in the following calculations to simplify the algebra, but we emphasize that SM^ may not 
necessarily be equal to 5M^. By differentiating Eq. (38), we obtain a general expression for the variation 
of the g'-factor 

^9nQm _ ^rrip / Cu 6M^ c^ 6Md \ 

5nqm rup \cu + Cd Mu Cu + Cd M^ J ' 

The proton mass, rup, and Mu^d are functions of the fundamental constants, and they can be formally 
written as nip = 'mp{vi,V2, ■ ■ ■ , Vn) and Mu,d = Mu,d ("Wi , '"2 , • • • , Vn), where the f j's are fundamental 
constants including m^, m^, rris, Aqcd, etc.. Therefore, Eq. (39) becomes 



Sgi 



NQM 



5nqm 





Vi drup 


^.- 


nip dvi 


n „ 

>r^ Svi 

- } Ki. 


i=l 





Vi 3Mu Cd Vi dMd 



Cu + Cd Mu dvi Cu + Cd Md dvi 



(40) 



This is our key equation in studying the dependence of the g'-factors on fundamental constants in the 
NQM approach, and the problem amounts to finding the expressions for rrip (f i , • • • ,Vn) and Mu^d {vi,- " j'^r. 
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3.1.1 The current quark mass and Aqqd dependence of rrip 

To get the coefficients of -^, that is, the first term in the square bracket of Eq. (40), we follow the 
procedure of [51, 52, 53], by defining Bq {q = u, d, s) and the 7r-nucleon sigma term, S^, in terms of 
proton matrix elements, 

mqBq = {p\mqqq\p) = rriq-^—^, (41) 

— CjfTi 

^TTN = {p\rn{uu + dd)\p) = m—^ , (42) 

om 

where rh = ^(rriu + "T^d)- The latter equalities of the above two equations come from the Hellmann- 
Feynman theorem [54] as noted by Gasser [55]. 
By using the strangeness fraction of the proton, 

_ 2B, ^ -^ _ ^L ^43^ 

Bd + Ba Stm ' 

where gq is the shift in the nucleon mass due to nonzero quark masses, and a relation from the energy- 
momentum tensor trace anomaly [56] for the baryon-octet members [57, 51, 58], 

Bd - Bs ms+ + ms- - rup - nin 
we can derive from Eqs. (41) and (42) the current quark masses dependence of rup, denoted as /x-q's, as 



/Tu 



-P 1 + :^ 1 + 1^ 



/T. = ^^^ = , '"^r, , . , (45) 



JTs = —- — 



mp 1 + ^ 



where 

Bd _ 2 + y{z- 1) 

5u 2z - y{z - 1) 
Motivated by the trace anomaly expression for rrip, 



(46) 



rup = rrinBu + rudBd + rusBs + gluon term , (47) 

we will write the remaining fundamental constants dependence of nip as 

which is the coefficient of (5Aqcd/Aqcd in Snip/mp. The argument behind Eq. (48) is the following: 
the gluon term has its origin in the strong interaction, and Aqcd, which is approximately the scale at 
which the strong interaction running coupling constant diverges, is the only mass parameter of the strong 
interaction in the chiral limit ni^ = m^ = nis = 0, and therefore in this limit all of the other finite 
mass scales of the strong interaction phenomena, including pion decay constant, the spontaneous chiral 
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symmetry breaking scale, etc., are related to Aqcd by some pure number of order one [59]. Note that the 
heavy quark (c, b, t) masses do not explicitly appear in Eq. (47) as discussed in [57]. Then as the only 
other variable besides the light current quark masses, we get Eq. (48) for the Aqcd dependence in nip. 
We note that the /xq's and the fx^ are also needed in the next section when we vary the electron-to-proton 
mass ratio, jj, = m^/mp. 

In calculating the /t^ 's and the Jt^ , we take the central values given in [60] for the current quark 
mass ratios, ^^ = 0.553 and ^^ = 18.9, the central value of the vr-nucleon sigma term suggested 
in [53], S^ = 64MeV, and we take ctq = 36MeV [61, 62] and rup = 938.3 MeV. The results are 

/t, = 0.027, /Td= 0.039, /t, = 0.363, /^^ = 0.571. (49) 

In the isospin-symmetric limit such that ttIu = rrid = m, which will be needed in subsections 3.2 
and 3.3, Eqs. (45) and (48) take simpler forms, 

jTm - , m = -T, ) jTg = 1 - jTm - IT, ■ (50) 

mp Imp '^ 

In calculating the values for this isospin-symmetric limit case, we take ^ = 25 [63], and the results are 

/t™ = 0.068, /7.^= 0.373, /^^ = 0.559. (51) 

3.1.2 Expressions for Mu d without an explicit quark sea 

To get the coefficients of -jf^ and -^, we need to model the constituent quark masses. Intuitively, Mu,d 
can be written as 

Mq = mq + aq^intAqcD (? = u,d), (A) (52) 

where aq^int's are pure dimensionless numbers. The argument behind this form is the following: if the 
strong interaction were switched off, the constituent quark mass would be identical to its corresponding 
valence current quark mass which is obtained from the electroweak symmetry breaking. On the other 
hand, in the chiral limit, w-u = rrid = rus = 0, the strong interaction is responsible for the entire con- 
stituent quark mass. The above intuitive expression for the constituent quark masses does not explicitly 
take into account the sea quark contribution, which if included will depend on the current quark masses, 
similar to the terms mqBq in the proton mass trace anomaly formula Eq. (47) [51]. However, one could 
argue that the sea quark contribution is already included implicitly in the second term of Eq. (52) to- 
gether with the virtual gluons contribution, since the dynamics of the quark sea and virtual gluons are 
determined by strong interaction, which is characterized in the second term. 
From Eq. (52), we obtain the coefficients of -jf^ and -j^ as 

TOu 9Mu _ rwu nid dM^ _ mg dM^ _ Aqcd dM^ _ m^ 

Mu drrin ~ Mu ' Mu dm^ ~ M^ dm^ ~ ' Mu Mqcd ~ Mu ' 

md 9Md _ rwd mu t^Md _ m^ dM^ _ Aqcd _9M^ _ ma 

Md drriA ~ Md ' Md dm^ ~ Md 9ms ~ ' Md OAqcd ~ Md " 

In calculating the above coefficients, we will use ^ = 0.553, the central value of md = 9.3 MeV in the 
modified minimal subtraction (MS) scheme at a renormalization scale of 1 GeV [60], and we will choose 

Mu = Md = 335MeV. 

3.1.3 Expressions for Mud with an explicit quark sea - linear form 

A method explicitly taking into account the sea quark contribution can be traced back to the internal 
structure of the constituent quarks [50]. Then, for ?71p(Aqcd, ^n^m^i^di^is) a linear realization of this 
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method is 

Mn = aiinAqcD + feu,lin"T-u + hd.lin^T-d + &s,lin"T'S , 

Md = aiinAqcD + &d,lin"lu + feu,lin"^d + Km^s , (B) (54) 

where we have related the coefficients in Mj with those in M^ following [50]. The coefficients aun and 
6q lin's are pure numbers. Each of the four terms of Mu ^ can be obtained by inserting Eq. (54) into an 
expression for the NQM based proton mass, rrip nqm, which we will discuss shortly (see e.g., Eq. (57), 
(59) or (61)). Then, applying the Hellmann-Feynman theorem 

drriq dM^ dniq dM^ dniq '^ 

An example of the application of the Hellmann-Feynman theorem within the NQM is given in [64]. 
From Eqs. (54) and (55), the coefficients of 5vi/vi (Vi = ?Ttu^d,s) Aqcd) of (5Mu/Mu and dM^^/M^ can 
be obtained as 



md dM^ ku mp/r^ - fed mp/r, (^ 



(56) 



Mu 9md Mu {kl - kj) 

Mu 9ms ~ Mu (feu + fed) ' 
Aqcd gMu _ 1 _ V^ "^q dM^ 
Mn Mqcd ~ ^ ^u 9mq ' 

where feu,d = ^7mIT ' ^^^ 7^^ (^^ = "^".d,s, Aqcd) are obtained from the corresponding ^^ 
by switching Mu -H- Md and feu -H- fed. 

To get feu,d» we consider the following NQM based proton mass formulae as examples. To zeroth 
order, the proton mass is the sum of the masses of its two constituent u quarks and one constituent d 
quark 

mp = 2Mu + Md , (57) 

so that 

feu = 2, fed = 1. (58) 

We will use M^ = M^ = ^ nip in Eq. (56) when Eq. (57) is taken as the NQM based proton mass 
formula. 

Without some interaction between the constituent quarks, hadrons with the same constituent quark 
compositions would have a same mass, a phenomenon which is not observed in nature. To break the 
mass degeneracy, a spin-spin hyperfine term is introduced [65], and the resulting proton mass is 

mp = 2M„ + M. + ^'(jlj-^), (59) 

where A' is a constant usually determined to allow an optimal fit to the baryon octet and decuplet 
masses [49]. This spin-spin hyperfine term is commonly attributed to one-gluon exchange [65], or, 
in the chiral quark model [66], it is explained as the interaction between the constituent quarks mediated 
by pseudoscalar mesons [67]. Although interpreted with relating to different degrees of freedom (gluon 
or pseudoscalar mesons) [58], this term nevertheless has a strong interaction origin, and therefore we 



ku = 2 + A'{^^^j:j-^^^], fed = 1 + T:rT72- (60) 
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will write the parameter A' as ahypA^^p, with a^yp a pure dimensionless number. From this formula, 
we get 

We will use M^ = M^ = 363 MeV and A' = (298.05 MeV)^ in Eq. (56) when Eq. (59) is taken as the 
NQM based proton mass formula. Note that we have tuned A' a bit compared to the value given in [49] 
to allow an exact fit to the proton mass. 

Eq. (59) can be further refined by adding to it the kinetic term of the constituent quarks and a con- 
stituent quark mass independent term Mq, which represents the contributions of the confinement potential 
and the short-range color-electric interaction [68, 64] 



4M2 MuMd / V ^u 2Md 



mp = 2Mu + Md + A" ( — ^ - — ^ ) + 5' ( -^ + ^ ) + Mo . (61) 



From the physical meaning of these two new terms, it may be reasonable to write the constants B' and 
Mq as CkinAg^jj and acccA^qcD. respectively, since the internal dynamics of a baryon is dominated by the 
strong interaction and the confinement is a strong interaction phenomenon. The constant A" needs to be 
re-fit after introducing the two new terms, and we write it as A" = a'^ K^^^. The parameters a\\^, acce 
and a[ are pure dimensionless numbers. We find that k^^^ from this formula are 

k-2 + A"(-^ l_\_^ feH-l + ^^--^ (62) 

We will use M^ = M^ = 335 MeV, A" = (4h76.4MeV)3, B' = (175.2 MeV)^ and Mq = 
—57.4 MeV in Eq. (56) when Eq. (61) is taken as the NQM based proton mass formula. Note that 
we have tuned Mq a bit compared to the value given in [64] to allow an exact fit to the proton mass. 

3.1.4 Mu d expressions with an explicit quark sea - NJL model 

As can be seen from Eq. (56), the explicit inclusion of the sea quark contribution in the linear form 
Eq. (54) encodes the information of both the NQM based proton mass formula and the /Tq's. However, 
different realizations from the NQM alone are also possible. Moreover, although the constituent quark 
masses are usually taken as fitting parameters in the study of hadron properties, it is certainly more 
illuminating if some concrete physical origin of these quantities can be given and encoded in their mass 
formulae. As suggested in [66, 69], the constituent quark masses are closely related to spontaneous 
chiral symmetry breaking. An example of the constituent quark mass formulae applying this idea is 
given by the three flavor Nambu-Jona-Lasinio (NJL) model [70, 71], where the constituent quark masses 
are obtained from a set of gap equations 

Mu = w-u - 2gs {uu) - 2go (dd) (ss) , 

Md = nid- 2gs (dd) - 2go {uu) {ss) , (63) 

ms - 2gs {ss) - 2go {uu) {dd) , (C) 



Ms = 



where {uu), {dd) and {ss) are the quark condensates which are the order parameters of the spontaneous 
chiral symmetry breaking, and they are calculated by one loop integral 

d'^v 1 

{uu) = -iNcTt ' ^ 



(27r)4 p-M^ + ie 
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and (^dd) i{ss)) is obtained by changing M^ to Mj (Mg). The Nc is the number of colors, and we take it 
to be the real-world value 3. This integration can be performed by introducing a three-momentum cutoff 
A3, and the result is 



3 

(uu) = -— ^Mu 



M^^[7Ws-M^J^i±^g±^ 



(64) 



The gs and g^ in Eq. (63) are the coupling constants of the effective four-point and six-point interactions 
of the quark fields in the NJL Lagrangian, and they are fixed, together with rrig and the cutoff A3, by the 
meson properties as explained in [64, 71]. We simply quote the result given in [64] ^ 

ms = 135.7MeV, c/^AI = 3.65, ^^A^ = -9.47, A3 = 631.4 MeV, (65) 

which we will use for our calculation. The other parameters we need in order to solve Eq. (63) are the u 
and d current quark masses, which we take m^ = ma = 5.5 MeV following [64]. Note that the form of 
Eq. (63) requires rriu = rrid if we assume M^ = Mj. The cutoff A3 characterizes the spontaneous chiral 
symmetry breaking scale, while the latter is related to Aqcd, as we explained in the paragraph below 
Eq. (48). Therefore, we will write A3 = ac,NjLAqcD, gs = as,NJLAQCD and g^ = ac,NjLA~^D. where the 
coefficients are pure dimensionless numbers. With these inputs, the constituent quark masses are solved 
from Eq. (63), with the values M^ = M^ = 335 MeV and Ms = 527 MeV, and we get 

mu 9Mu _ rwd dM^ _ m^ dM^ _ m^ dM^ _ 

Mu drriu Mj dm^ ' ' M^ dm^ M^ drria 

iris dMa iris QMa Aqcd dM^ Aqcd QMa 

J ^ = TT^^r^ = 0.0628, ?T ^^ = ?T ^^ = 0.8947. (66) 

Mu drris M^ orris M^ 9Aqcd M^ SAqcd 

3.1.5 Results and discussion of the NQM approach 

We can now calculate the dependence of the nucleon magnetic moments on mu,d,s and Aqcd, by 

O^NQM _ 5nQM O^NQM (f^l\ 

9cxp dexp 5nqm 

where g^xp is the measured value of the (/-factor, which equals 5.586 for proton, and —3.826 for neu- 
tron [63]. The first term in the square bracket of ^^'^^^ (Eq. (40)) is given in section 3.1.1, while 
the second term in that square bracket can be obtained from section 3.1.2, 3.1.3 and 3.1.4 for each of 
the three different constituent quark mass models we have considered. The calculated coefficients of 
-^ {vi = mu^d,s) Aqcd) of ^'^'^^'' for proton and neutron are listed in Table 3, where the constituent 
quark mass formula used for each row is labeled as A, B or C, representing Eqs. (52), (54) or (63), 
respectively, while the 1, 2, or 3 following the label B represents Eq. (57), (59) or (61), respectively. 
In all of the cases listed in Table 3, we use the /t given in Eq. (49). Note that for case C, there is a 
slight inconsistency due to our choice of m^ = m^, though this has only a minor numerical effect on the 
resulting k's. 

The coefficients in Table 3 show a relatively strong dependence on the constituent quark mass models 
used. Most of the coefficients in A and C are closer and much larger compared to their corresponding 
values in B. While the first and the second terms in Eq. (40) are independent of each other for A and C, 
the same /^q's appear in both terms of Eq. (40) for B, as can be seen from Eq. (56) and thus these two 
terms are largely canceled due to a relative sign. We can also see a relatively strong dependence of the 
coefficients on the NQM based proton mass formulae when comparing the rows Bl, B2 and B3. 



'We have tuned the values of Qs and go relative to the values given in [64] to allow Mu — Md = 335 MeV and A4s 
527 MeV as exact solutions of Eq. (63). 
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Table 3: The coefficients kj of -^ {vi 
(right); see Eq. (40) for their definition. 
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Figure 1 : The dependence of the coefficients in A for the proton on S,, 
B3 (right). 



(left) and of the coefficients in 



Furthermore, there is an uncertainty in the coefficients listed in Table 3 due to the uncertainty of the 
vr-nucleon sigma term S^. A discussion of the impact of the uncertainty of Tjj^^ on the interpretations of 
experimental searches for dark matter can be found in [53]. We plot the dependence of the coefficients in 
row A for the proton on I],n^ in the left panel of Fig. 1 . A similar plot of the coefficients for the proton in 
row B3 is given in the right panel of Fig. 1. As can be seen from these plots, the coefficients of Snis/nis 
and 5Aqcd/Aqcd show a strong dependence on the value of S^. Therefore it is important to pin down 
the value of S^ if this quantity is used in the study of the current quark mass and Aqcd dependence of 
the proton (7-factor. The same conclusion applies for the neutron gr-factor, for which the behavior of the 
plots are similar to that shown in Fig. 1. 

In addition to the relatively strong dependence of the kj on the proton and constituent quark mass 
formula as well as the value of S^ we have discussed above, some other comments for this NQM ap- 
proach in the study of the dependence of the nucleon (7-factors on the fundamental constants are in order. 
Our assumption that the various parameters in the constituent quark mass formulae and the NQM based 
proton mass formulae take power law forms for Aqcd may be valid only in the chiral limit. Therefore, 
some current quark mass dependence may be lost and the Aqcd dependence may not be very accu- 
rately determined from these formulae. To get a more accurate dependence, one may also wish to con- 
sider relativistic corrections [72] and/or corrections based on higher-dimension terms in the chiral quark 
model [66] for Eq. (37), and then the dependence on the current quark masses and Aqcd will change 
correspondingly. Finally, in the above analysis, we did not consider the electromagnetic contribution to 
the proton mass or the constituent quark mass formulae, and thus we may have missed some dependence 
on the fine structure constant in this approach. 
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3.2 The chiral perturbation theory approach 

The second approach we consider is xPT, which provides a systematic method of addressing the low 
energy properties of the hadrons [73, 74]. In contrast to the strong model-dependent NQM approach 
we considered in the previous subsection, xPT can give model-independent calculations of the nucleon 
magnetic moments within a perturbative field theory framework in terms of the hadronic degrees of 
freedom. However, as we will see, our goal of extracting the current quark mass and Aqcd dependence 
is limited by our lack of knowledge of the accurate values of the coupling constants, the so called low 
energy constants (LECs), appearing in the effective Lagrangians of xPT. These Lagrangians, and the 
Feynman diagrams generated by them, are organized according to a power counting scheme, and the 
number of LECs we will have to deal with increases as we include higher order contributions to the 
nucleon magnetic moments. 

By construction, the LECs in the SU(3) xPT which we will consider do not depend on the light quark 
(u, d and s) masses, and they should in principle be calculable in terms of the heavy quark (c, b and t) 
masses and Aqcd- Without the ability to solve non-perturbative QCD, the LECs are usually determined 
by fitting to experimental data for the pertinent physical observables, or estimated theoretically by QCD- 
inspired models and some other approaches (e.g., the resonance saturation method), and they can also 
be fixed by lattice calculations (for a discussion of the LECs, see for example, [75] and the references 
therein). Most of the LECs are renormalization scale dependent in such a way that they cancel the 
renormalization scale dependent loop integrals so that the final results for the physical observables are 
renormalization scale independent. Furthermore, the values of the LECs are expected to be given by 
dimensional analysis [66, 76] up to numerical factors of order one. Since the two quantities involved 
in such an analysis, namely, the Goldstone boson decay constant (for the meson octet) and the typical 
mass of the light but non-Goldstone states, are both pure numbers times Aqcd in the chiral limit, we will 
assume that all the LECs under discussion are functions of Aqcd and the renormalization scale, and by 
this assumption we neglect the heavy quark mass dependence in the LECs. 

For xPT in the meson-baryon sector, needed for the calculations of nucleon magnetic moments, 
there exist several renormalization schemes in the literature to ensure consistent power counting which is 
troubled by the introduction of the baryon mass as a new scale which is non- vanishing in the chiral limit. 
Among these renormalization schemes, the most studied in the early days in the calculations of octet 
baryon magnetic moments is the heavy baryon chiral perturbation theory (HB^PT) approach [77]. Due to 
a strong cancellation between the leading order 0{q^) and the next-to-leading order 0{q^) results for this 
approach {q denotes external momentum in the power counting scheme), one is forced to consider still 
higher order contributions. We will consider the results for this approach to order 0{q'^) [78, 79, 80, 81], 
with (HBwD) and without (HBw/oD) the explicit inclusion of the baryon decuplet states in loops. We will 
also consider a result from a more recently developed extended-on-mass-shell (EOMS) renormalization 
scheme [82], which gives more convergent results at 0{q^) without [83] or with [84] the inclusion of 
decuplet states in loops. We will restrict our attention to the EOMS without decuplets to avoid the 
introduction of several new parameters which do not improve the convergence. 

At leading order, the octet baryon magnetic moments can be calculated from the Feynman diagrams 
of chiral order 0{q^), and the results for both the HB^PT and EOMS approaches have the same expres- 
sions as linear combinations of two LECs ii^ and ^p, 

I^B = «B = "B/^D + "S^F , (68) 

where a^ = 1/3 and a^ = 1 for the proton, and a^ = —2/3 and q^ = for the neutron. Note that 
we are writing down the magnetic moments (rather than the anomalous magnetic moments) directly in 
units of fi^, and therefore the jj^p value we use may differ by 1 compared to the value given in some of 
the references. 

At 0{q^) and higher order, the results of HB^PT and EOMS differ. For the HBw/oD approach, we 
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use the result given in [78, 79, 80, 81], essentially using the notation of [80]. At 0{q^) it is 



where 



2 

'3* 



I3; = -{D + Ff, f]^ = --{D' + 3F^), Pl = {D + F)\ ^^ = -{D - Ff 



D and F are dimensionless LECs. We will use the empirical values nip = 938.3 MeV, the pion decay 
constant /^r = 93 MeV and the kaon decay constant Jk = 1.2/7^ in our calculations, but we will consider 
these quantities take the forms of their corresponding LECs, namely, the average octet baryon mass, ttiq 
and the Goldstone boson decay constant Fq, since the differences between these quantities and their 
corresponding LECs give contributions to the octet baryon magnetic moments beyond 0{q^) which 
is the highest order we will consider. Therefore, in contrast to the NQM approach we considered in 
the previous subsection, we do not need to vary the proton mass in this subsection in calculating the 
current quark masses and Aqcd dependence of the nucleon magnetic moments. For the same reason, 
we take the pion mass, M^, and kaon mass, Mk, at their empirical values of 138 MeV and 495 MeV, 
respectively, in our calculations, while we only take their lowest order forms M^r = {23^1%)'^/'^ and 
Mk = [BQ{rh + rris)]^''^, where the Bq is an LEC with mass dimension, when considering light current 
quark masses and Aqcd dependences. We will work in the isospin-symmetric limit in this and the next 
subsection such that m^ = rud = m. The expression for the 77 mass, M,, = [2/3Bq (m + 2?7is)] ^' ^ is also 
needed at 0{q*) in the HB^PT approach and at 0{q'^) in the EOMS approach, and we use its empirical 
value of 548 MeV. 

For the case HBwD, we have the following terms in addition to Eq. (69) [78, 80, 81], 

^'B= E ^i^(Mx,<5,A)/3f , (70) 

X=-K,K ^•'^ 

where 

2 f 2VM2 - 52 [7r/2 - arctan {5/VM'^ - 6'^)], M -^ 5, 

7rF{M,5,X) = -6ln—^ + 
A^ 



and 



-2V52-M2ln[((5 + V<52-M2)/M], M < 6, 



o/TT _ __r>2 niK _ _r>2 ol-n _ ^nl oiK _ \p2 



where C is a dimensionless LEC, A is the renormalization scale in dimensional regularization, and b is 
the decuplet-octet mass splitting for which we take to be a number times Aqcd with a value of 300 MeV. 
At 0{(f) in the HBxPT approach, more LECs appear in the results and the formulae become lengthy. 
For the case HBw/oD, we take [79] ^ 

(4) (4,c) , (4,d+e+/) , (4,9) , (4,/H-J) , (4,j) ,„,. 

f^B=f^B +f^B +IJ>B +IJ'B +fJ'B^ (71) 



with 



(A r) 111 ^4 \ 2 2 1 

l^p ' =03 + 04 + -05 + -ae - -aj, < ' ^ = --as - -qq - -07, 



^There is a misprint in the third term of fi^ ' "^ in [79], where the sign in front should be '- 
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M2 



(4,d+e+f) 






X = TT,K 



X 



.ln^+ E 



A 



Vb 



X Ml 



Y^ VB 

X=n,K,r) 



X Ml 



SvrVi 



X=TT,K,-q 



, Mx 

31n— ^ + 1 



levrVi 



as, 



m— h 1 

A 



17 



with 



and 



with 






-^D — MF, 



cif 



/^D — ^^F, 



1 



7^ = -(Z) + F)^(^B-^^) 



^1 



K 



vl 



e 



P)I1 



--p-3F)2(^^ + 3^i.), 



^D2 - 2Z)F + F2)^^ - (Z) - F)2^f , 
y 



Vn 



1d^ + ^DF + f2)^z) + (I? - F)Vf, 



< 



^(Z)-3F)Vd, 



-p + i^)^ 



6^ 

p,n 



.g^^^ 



-Z)^ 



df + -f'^, 

2 ' 



p,n 



1 

12 



(D - 3F)2, 






X=7r,_ft: 



47r2/' 



X 



A 



2a9 + 2 aio + 



8m 



i: 



-2ag - 2 aio + 



Srrir 



K 



7f 



+ 4 aio + 



Sniy 



-2a9 + 2 oio + 



Srrir 



and 



and 



with 



(4,h+j) 
Mi? 



Mg"' 



X=7r,X 

- E ^ 

X=-K,K 



Ml 
16vr2/2. 



X 



rrir 



^ 27r2/l 



51n^ + l 



21n^ + l 



^^ = (D + Ff [Mlaii + {Ml - Ml) ai2] , ^f = - [(3F + Of M^ + 3(Z) - FfMl 
ei = -{D + Ff [Mlau + {Ml - Ml)ai2] , 0^ = {D - FfMlan. 



an, 



In the above formulae, the LECs 08,9,10,11,12 (labeled 69,10,11, d,f in [79]) with their values in units of 
GeV~^ are as = 0.81, ag = 0.95, aio = 0.36, an = -0.192 and ai2 = -0.210, where the first three 
are estimated by the resonance saturation method which takes into account the contribution from the 
baryon decuplet while the other two are determined by fitting to the baryon octet masses, as explained in 
detail in [79]. We take the value of the rj decay constant to be /^ = 1.2/,^, but we will consider it taking 
the form of its corresponding LEC, Fo, for the same reason explained above for the other Goldstone 
boson decay constants. 
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At this order, an uncertainty arises for our extraction of the dependence of magnetic moments on 
the current quark masses and Aqcd due to the uncertainties of the values of the LECs appearing in the 
counter term Lagrangian, and they are denoted as 03^4^5^3,7 in Eq. (71). These dimensionless numbers 
have actually absorbed the light current quark masses in the counter term Lagrangian in contrast to their 
corresponding true LECs which are independent of the light current quark masses. Therefore these 
redefined LECs should contain a factor ms/AqcD. if we neglect the contributions from m as its value is 
much smaller than mg. Two other LECs in the counter term Lagrangian are also present at this order, 
and they are combined with the two LECs appearing in the 0{q'^) result Eq. (68), as 

IJ-D,F^ fJ'D,F + 4:Bo{2m + ms)flD,F, (72) 

where fj,D,F are LECs appearing in 0{q'^) counter term Lagrangian. Then all seven of these redefined 
LECs, /Ud,f and 03^4^5^6,7. are used as fitting parameters to perform an exact fit to the seven available 
octet baryon magnetic moments. Since they are used as fitting parameters, and indeed different values 
for them are obtained with and without the explicit inclusion of baryon decuplet states in loops, as well 
as when different values of other LECs are used for the fittings (see the discussion below), it is hard to get 
an accurate extraction of the light current quark mass and Aqcd dependence from these redefined LECs. 
For the light current quark masses dependence, we will only consider the rris dependence for a3^4^5^6,7» 
while we will not try to extract such dependence for the redefined ij.d,f (denoted as ai^2 in [80]), since 
we do not know the relative size of the two terms on the right hand side of Eq. (72), where only one of the 
two terms has the light current quark mass dependence, although such dependence in these two redefined 
LECs may be not small, as suggested in [79] when comparing the fitting values up to 0{q^) with the 
ones up to 0{q'^). We take the values fi£) = 3.71, hf = 3.25, 03 = —0.50, 04 = —0.15, 05 = —0.25, 
oe = 0.50 and aj = —0.46 given in [79], where F = 0.5, D = 0.75 is used, and the renormalization 
scale A is taken to be 0.8 GeV. 

For HBwD at ©(g^), we take [80] ^ 



X=n,K,rj ■' ^ 



+ E 



1 



^— ^ 327r2/2 

X=-K,K,r) •' ^ 



(7^^ - 2 ,^'# aB)L(3/2) (Mx , <5, A) + 7^^ L'(3/2) (Mx , <5, A) 



(73) 



where 



and 



with 



M2 

L(3/2) (M, 5, A) = m2 In^ + 2ti5 F{M, 6, A), 



L'(3/2)(M,<5, A) = M^ln^ + ^G(M,<5, A), 



M2 

ttG{M, S, A) = -5^ In— ^ + ttM^ + 
A^ 



-2(M2 - (52)3/2 [^/2 - arctan {S/VM'^ - 6'^)], M ^ S, 
-2{6^ - M2)3/2 in[(5 + V(52 - M2)/M], M < S, 



'There is a misprint in the form of the -kG(M, S, A) for the case M ^ S, where the sign in front should be ' — '. 
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and 



''-'■- 7f = ^CVc, %^ 



l'; = ^c'^lc, %'' = ^c'fic, 7';> = o, 



7n 



27 



-c'fic, i':^ = --c'nc: i:? = o, 



7, 



P 9 

In 






'-- ^C{D + F)^1T, j'f = -^CFfiT, ln^ = 0, 



Pill 



2C2, 



Pill 



c^ 



^^'^n = 0. 



The other coefficients are related to the ones given in Eq. (71), as 

^B +Vb ' 7b 



lTT,K 

7b 



V 



IX 
'B 



3(pB, 



,(4,c) 



and fJ-)^'"'' is the same as the case of HBw/oD. 

In this case, since different LECs are used as inputs for the fittings in comparison to the HBw/oD 
case, the resulting fit values of the seven redefined LECs are different, and we take their values from Case 
(b) in Table II of [80], oi = 3.946, 02 = 2.353, 03 = -0.001, 04 = -0.172, 05 = 0.569, oe = 0.694 
and 07 = —1.165, corresponding to the LECs inputs F = 0.5, D = 0.75, C = —1.5, fir = —7.7 and 
He = 1.94. A renormalization scale A = 1 GeV is used. 

For the EOMS approach, to minimize the number of LECs involved and thus perhaps the uncertain- 
ties introduced by them, we only consider the result given in Eq. (2) to Eq. (5) of [83] which does not 
include the baryon decuplet states in loops. The result is up to 0{q'''), and the loop integrals are finite 
The values ^u/) = 3.82 and fip 
are determined by a fit to minimize the x^ = Yl (/"th 
quantities in the formula, F^ = 1.17 f^^ is the average of the physical values of /t^, /k and f^j, and we 
still use Mb = 938.3 MeV, /t^ = 93MeV, and the same values for M^^^k,?; as specified above. We take 
D = 0.80 and F = 0.46 as used in [83] for this EOMS approach. 

We list the results of the two HBxPT and the one EOMS approaches in Table 4, where we also need 
to specify the ratio of rris to m, for which we use 25. Note that as we discussed above, we have assumed 
that the LECs (ai^2,3,4,5,6,7» correspond to the original LECs before the re-definition) depend only on 
Aqcd and the renormalization scale. A, i.e., they do not depend on the light quark masses and there is 
no dependence on the renormalization scale in the full result. Therefore, the coefficient of (5Aqcd/Aqcd 
must be equal and opposite to the sum of the light quark mass contributions. 



2.20 denoted as b^ and b^ in [83], after performing the EOMS scheme, 
/f^cxp)^ as explained in that reference. For other 



Table 4: The coefficients, kj, of 
(left) and the neutron (right). 



5vi 



(Vi 



m, m^ 



A 



QCD. 



of -^2s£T ^ defined as in Eq. (40), for the proton 

ycxp 



ZitiiYii — '^'^d 



I^QCD 



Ztrxjyi^ — Zirvif^ 



'^QCD 



HBw/oD 

HBwD 

EOMS 



-0.050 

0.034 

-0.049 



-0.50 

0.17 

-0.031 



0.54 
-0.21 
0.080 



-0.16 -0.14 0.30 

-0.050 0.32 -0.27 

-0.11 0.014 0.097 



We see from Table 4, the numbers in each column differ considerably for nucleon magnetic moment 
formulae from different renormalization schemes and depend on the explicit inclusion of baryon decuplet 
states in loops. As we mentioned in the beginning of this subsection, we believe this discrepancy comes 
in a large part from our lack of knowledge of the accurate values of the LECs. In particular, many of the 
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LECs involved in our calculations are used as fit parameters for the octet magnetic moments and masses, 
while their true values may be quite different from the numbers obtained by these fits. For example, if 
we use the values for Case (a) F = 0.4, D = 0.61 and C = — 1.2 in [80], the resulting values of the 
a's also give an exact fit to the seven available octet baryon magnetic moments, with a prediction for 
the EA transition moment similar to the one given by Case (b). However, one can see that many of the 
corresponding a's for Case (a) and Case (b) differ greatly, and indeed, for Case (a) we get the coefficients 
from left to right of Table 4 as 0.014, 0.14 and -0.16 for the proton, and -0.039, 0.26 and -0.23 for the 
neutron, which are different from the results of Case (b). Therefore, it is crucial to pin down the values 
of LECs before one can make a better extraction of the light quark masses and Aqqd dependence in the 
XPT approach. 

One can also estimate the dependence of rup on the current quark masses and Aqcd from a formula 
for nip within xPT. Such dependences can be used when one varies the electron-to-proton mass ratio, 
fi = nic/mp. However, as we explained at the end of the previous subsection, we prefer to use a common 
set of values for the coefficients of Hq and kqcd of Snip /nip. Those values for the isospin-symmetric 
limit case are listed at the end of section 3.1.1. 

3.3 The approach combining xPT and lattice QCD 

As another approach to study hadronic physics, lattice QCD provides a promising way to extract the 
current quark masses dependence of the nucleon magnetic moments, because one can do explicit cal- 
culations by assuming a sequence of different current quark masses in lattice computations, although in 
practice the computational cost is a limitation. Since most of the current lattice computations are still 
using input current quark masses much larger than their empirical values, an extrapolation of the lattice 
results to the physical point is needed. In the extrapolations for the physical observables, terms hav- 

1/2 

ing non-analytic behaviors, rriq and ni^ log rriq, etc., which are predicted by ^PT and have important 
contributions near the chiral limit, must be considered. 

An earlier study of this combined lattice and ^PT approach for the nucleon magnetic moments uses 
an encapsulating form which is the Pade approximant [85], 

^ JM^) = ^ — '^ , (74) 

Mo ^ 

where Xp,n are fixed by the leading non-analytic term given by xPT, while /Uq and c are allowed to vary 
to best fit the lattice data. 

A later development takes the finite range regulator (ERR) [86] as the regularization method rather 
than the traditional dimensional regularization for the results we discussed in the previous subsection, 
and the cut-off parameter in the ERR is a mass scale which can be interpreted as the inverse of the size 
of the nucleon. 

The current quark masses dependence for the nucleon magnetic moments is given in [28], and we 
simply quote the result there without going into any detail 

(Jon Srh 6ni^ (Jo,, Sm 5m^ 

JlE = -0.087 0.013 -, -i^ = -0.118 +0.0013 ^. (75) 

Qp m iris gn m m^ 

As the same argument we made for the xPT approach in the previous subsection, all parameters without 
light quark masses dependence are either pure numbers or are pure numbers time Aqcd- Therefore, we 
obtain 

^ = -0.087^-0.013^ + 0.100^^^, 
5p m nis Aqcd 

6q„ 5m Snis 5Aacr> 

-^ = -0.118 +0.0013 +0.1167 , ^ . (76) 

gn m nis Aqcd 
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4 Atomic clock constraints 

4.1 Methodology 

As we have seen in section 2.4, the frequency shift is related to {gp, g^, b, fi, a} by the relation 

^ = V^ + A,„^ + a4 + a/ + A.«, (77) 

where the coefficients A are given explicitly in Table 2. Then, in section 3, we have expressed the 
dependence of the g'-factors as 

6gp Sruu , 5md , 6ms (JAqcd ,„„, 

5p "T-u rud ^ rris ^ Aqcd 

6gn Sniu ^ 5md , 5ms <JAqcd ,„„, 

5n "T-u rrid ms Aqcd 

where the coefficients Kj have been calculated for different models and collected in Tables 3 and 4 and 
Eq. (76), as well as the dependence of the proton mass 

^ = /T.^ + hj-^ +fTf-^ + hj-^, (80) 

mp mu md m^ Aqcd 

where the /t^ are given in Eqs. (49) and (51). Also, following [27, 31], b depends on the quark mass and 
Aqcd, and there it is found 

Sb 5m 5 Aqcd .„., 

y=7.- + 7ccD^^^, (81) 

with 

7q = -7qcd = -0.11. (82) 

Assuming for simplicity that all Yukawa couplings are varying similarly, i.e., 5hi/hi = 5h/h, the 
expansions (78), (79), (80) and (81) can be inserted in Eq. (77) to obtain 

''^^ - X ^ ^ X "^ ^ X ^QCD t a 

— A/j— + Ay h Aqcdi r Aq — . {pi) 

VAB h V Aqcd a 

The coefficients A are easily computed to be given by 

•^fe = V('^Up + '^dp + Ksp) + XgAl^Un + /^dn + Ksn) + Ab7q + A^(l - /t, - /Td - /Ts)(84) 

At, = A/i (85) 

Aqcd = AgpKQCDp + Ajj„KQCDn + Afe7QCD - A^/Tg (86) 

Aa = Aq. (87) 

The form (83) makes no assumption on unification and only relies on the fact that all Yukawa couplings 
are varying in a similar way. It is important to note here that the dimensionality constraint on the fx^ , kj 
and 7j parameters implies that Aqcd = — A^, so that Eq. (83) actually depends only on the combination 
of X = /iv/Aqcd and a as 

^ = a4 + A„«. (88) 

This would not be the case if we had not assumed that 5hi/hi = 5h/h for all Yukawa couplings. 
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Our first hypothesis concerning unification allows one to express the variation of the QCD scale by 
means of Eq. (2) so that 

^h + TT'^QCD ] T ^ { ^v ^ TT'^QCD ^ ( Aq, + RXqcd 



J h \ 9 J V \ J a 

( h v\ a 

= Hs\- + -\ +H^-. (89) 




The second hypothesis on unification assumes that the variation of v and h are related by Eq. (4) so that 

Aft + -Aqcd ] {^ + S)— + [Xa + RXqcu 



/ n \ /a 

= Hs{l + S)^ + H^-. (90) 

n a 

The last hypothesis assumes that the variations of h and a are related by Eq. (5) so that 



l^AB 



2\^h+ g Aqcd ) (1 + 5) + ( A« + i?AQCD 



a 
a 



= {\hs{1 + S) + H^)-^Co.{R,S)-. (91) 

The two last equations define the parameter Ca{R, S). 

The forms (89-91) imply increasing assumptions on the unification mechanisms and are thus be- 
coming more and more model-dependent with the advantage of reducing the number of fundamental 
constants, hence allowing one to draw sharper constraints from the same experimental data. 

The coefficients introduced above can be easily calculated from Table 2 for the coefficients A, Ta- 
bles 3 or 4 or Eq. (76) for the coefficients Ki, Eq. (49) or Eq. (51) for the coefficients fx^, and Eq. (82) 
for the coefficients 7j. As an example, we provide the value of the coefficients Ca assuming S* = 160 
and R = 30 for the 9 models considered in this article. It is important to stress that this coefficient is 
almost always larger than one and typically of order 5 - 30 in absolute value. 

We can check that the effect of varying the nuclear radius is indeed much smaller than varying the 
other parameters. This effect can be included by adding a term er{rh/m — Aqcd /Aqcd) to Eq. (77). 
Using the values listed in Table IV of [31], we have e^ = —0.004 for the Cs-Rb clock system, while 
er = —0.007 for the other five clock systems involving Cs. These amount to an adjustment of —0.13 in 
the numbers in the first column of Table 5, and —0.23 in the other five columns. 

4.2 Single Experiment Constraints 

We are now in a position to combine our results for the dependence of the nucleon (/-factor on funda- 
mental parameters with the limits imposed from atomic clock measurements. For each experiment, we 
can derive a limit on the variation of the fine structure constant under a number of sets of assumptions. 

For example, assuming first that the only dependence of uab on a is related to the coefficient X^ 
(i.e., we assume that gp, g^, b and /x remain constant), we can use Table 2 to extract a limit on a/ a for 
each experiment from 

^ = ^^. (92) 

a Xa ^AB 

In contrast, when we take into account the contributions from coupled variations, and we assume the 
relation between vab/^ab and a /a given by Eq. (91) we obtain simply 

- = yr^^- (93) 

a L-a VAB 
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Table 5: The coefficient Ca assuming S" = 160 and R = 30 for each of the models for the nucleon 
magnetic moment and for the various combinations of clocks discussed in this article. 





Cs-Rb 


H-Cs 


Hg-Cs 


Yb-Cs 


Sr-Cs 


SFg-Cs 


A 


-16.53 


13.86 


17.06 


12.96 


13.80 


4.56 


Bl 


-2.26 


20.16 


23.36 


19.26 


20.10 


10.85 


B2 


-6.79 


18.16 


21.36 


17.26 


18.10 


8.85 


B3 


-5.29 


18.82 


22.02 


17.92 


18.76 


9.51 


C 


-13.37 


15.26 


18.46 


14.36 


15.20 


5.95 


HBw/oD 


19.27 


29.33 


32.53 


28.43 


29.27 


20.22 


HBwD 


-8.57 


17.01 


20.21 


16.11 


16.95 


7.89 


EOMS 


0.49 


20.97 


24.17 


20.07 


20.91 


11.86 


XPT+QCD 


1.20 


21.29 


24.49 


20.39 


21.23 


12.18 



Thus, the improvement in the limit from each individual experiment due to the theoretical assumption of 
coupled variations is given by Cq/Aq. These factors are tabulated in Table 6 for each experiment and 
model for (7p n. 

Table 6: The enhancement factor Ca/Xa assuming S = 160 and i? = 30 for each of the models for the 
nucleon magnetic moment and for the various combinations of clocks discussed in this article. 





Cs-Rb 


H-Cs 


Hg-Cs 


Yb-Cs 


Sr-Cs 


SFe-Cs 


A 


-33.73 


4.90 


2.83 


6.72 


4.98 


1.61 


Bl 


-4.61 


7.12 


3.87 


9.98 


7.26 


3.83 


B2 


-13.86 


6.42 


3.54 


8.94 


6.53 


3.13 


B3 


-10.80 


6.65 


3.65 


9.28 


6.77 


3.36 


C 


-27.28 


5.39 


3.06 


7.44 


5.49 


2.10 


HBw/oD 


39.32 


10.36 


5.39 


14.73 


10.57 


7.14 


HBwD 


-17.48 


6.01 


3.35 


8.34 


6.12 


2.79 


EOMS 


1.00 


7.41 


4.01 


10.40 


7.55 


4.19 


XPT+QCD 


2.45 


7.52 


4.06 


10.56 


7.66 


4.30 



As one can see, there is a strong model-dependence on the resulting limits on a/ a. Overall the 
enhancements range from ~ 1 to ~ 10. For example, let us consider the case of the Cs-Rb atomic clock 
system. Ignoring the variations in all other constants, this clock would yield a result 



a 



-i-16„„-l 



- = (1.02 ± 10.82) X 10^'V" 
a 



(94) 



In contrast, coupled variations, according to the factors in Table 6, improve this result by as much as a 
factor of 39.32 using the HBw/oD model for ^fp „, yielding 



a 



16, ,^-1 



- = (0.03 ±0.28) X lO^^^yr 
a 



(95) 



Cases A and C also make substantial improvements in the limit for the Cs-Rb clock system. On the other 
hand, there is no gain for case EOMS, or even a weaker limit if the nuclear radius effect is taken into 
account. 
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4.3 Combined Experimental Constraints 

While the results of individual experiments can be substantially improved by coupled variations, two 
clock systems (Dy and Hg-Al) are independent of any assumption on unification and lead to model- 
independent limits on a. We next combine the available results to obtain a single limit on a for each 
choice of model for gp^^- 

Each of the eight experimental results used in this article can be written as 

— In UAB = VAB ± Sab, (96) 

at 

listed in Table 2. From a theoretical point of view, the expression for uab depends on a set of con- 
stants, X, chosen as being independent and on our hypothesis on unification schemes. If we assume 
dlnvAB{^)/dt — rjAB to be Gaussian distributed and all the experiments to be uncorrected, then the 
best-fit for the set of constants x is obtained by maximizing the likelihood, or equivalently by minimizing 






2^ N, X""^ [!^AS 



(97) 



AB "^S 



The 68.27%, 95%, and 99% confidence level (i.e., la, ^^ 2a and ~ 3a) constraints are then obtained by 

Ax^ = (1, 3.84, 6.63) if dim(x) = 1 and Ax^ = (2.30, 5.99, 9.21) if dim(x) = 2. 

4.3.1 Constraints on the QED parameters 

Let us start by assuming that {gp, g^, b, fi, a} are independent parameters. One can use the Hg-Al clock 
to constrain the variation of a and then use the six clock combinations that depend on the five parameters 
to set a constraint on {gp, gj^, h, fi}. However, from Eq. (30), we note that the ratio of the coefficients of 
6gijgn and 6b/b is gn/{gn — dp + 1)» which is independent of the clock systems we are considering. 
Also, from Table 2, we note that the value of Xg^/Xg^ for the Cs-Rb clock is very close to that of the other 
five clock combinations. Therefore, for the purpose of constraining the QED parameters, gp, g^ and b are 
not independent, and we can only constrain their combination, namely, gcs- The combined constraint on 
gcs and fi is depicted on Fig. 2. Note that if a different method in the calculation of gr-factors of ^^Rb 
and ^•^^Cs, and/or other clock systems, are used, such that the ratio A^p : A^^^ : Af, is not the same for 
different clock combinations, then gp, g^ and b can be taken as independent parameters. 

As we know from our analysis, such a hypothesis is not correct since the variations are expected to 
be correlated but this shows the result one would have derived without any knowledge on QCD. 

4.3.2 Constraint on a 

As in the previous subsection, we can consider first the constraint obtained using the form (91) that 
depends on 6a/ a alone. Minimizing x^ for a single variable is equivalent to taking the weighted mean of 
rjAB/Xa with an uncertainty Sab/Xo- This result can be compared with that assuming coupled variations 
using the coefficients Ca, given in Table 5. In this case, the weighted mean replaces A^ with Ca- 
in order to determine the effect of coupled variations, we compare the constraints arising from the 
combination of the eight experiments to the one obtained from the combination of 6 clocks (that is 
neglecting the Dy and Hg-Al clocks). The results are presented in Table 7 and shall be compared to the 
same analysis assuming that only a is varying (i.e., keeping gp, g^, b and /i constant). We find 

- = -(2.14 ± 2.30) X lO"^'^ yr"^ (98) 

a 
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Figure 2: Constraints on the variation of parameters {(7csi /f^} assumed to be independent once the con- 
straint from the variation of a from the Hg-AI clock is taken into account. Solid, dashed and dotted 
contours correspond to 68.27%, 95% and 99% C.L. 



for the combination of the 8 experiments and 



a 
a 



(5.24 ±6.40) X 10"^^ yr"^ 



(99) 



for the combination of the 6 experiments. We also remind the reader that the Hg-Al experiment alone set 
the constraint 

" - (1.65 ± 2.46) X 10"^'^ yr"^ (100) 



a 



which shows that there is little gain in combining the 8 experiments compared to this experiment alone. 
When 5(p, g^^, b and fx are allowed to vary in the combination of the 6 clocks, there is a gain of a factor 
of order 4 so that the constraint obtained from the combination of these 6 clocks assuming unification 
becomes as strong as the constraint obtained from Hg-Al alone. When combining the 8 experiments, the 
gain is less than a factor of 2, due to the fact that the limit arises mostly from the Hg-Al experiment which 
does not depend on gp n- These results are summarized in Table 7 and each result can be compared to 
the single Hg-Al result given in Eq. (100). 



4.3.3 Constraint on hv 

As a second application, we can use the constraint (22) arising from the Hg-Al clock to obtain a bound 
on the time variation of a that is independent of the other constants and then use the 6 other clocks to set 
a constraint on the combination of parameters hv, assuming the form (89) to set a constraint on Shv/hv 
alone. This requires the knowledge of the coefficients Hs and Ha and we assume that R = 30, but it 
does not depend on the coefficient S. 



The constraints for each model are summarized on Table 8. It ranges between 



(hv)- 



hv 



< 20.43 X 



10 



-16 yj. 1 ^^^ 



(hv)- 



hv 



< 13.52 X 10 ^^ yr ^, respectively for models A and HBw/oD and it turns out 



that the model-dependence for this constraint is mild. 
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Table 7: Constraints on the variation of a assuming the unification relation (91) and the values of Ca for 
S = 160 and R = 30. We compare the constraints obtained from the combination of the 8 clocks and 
the constraints obtained from the 6 clocks (i.e. without Dy and Hg-Al). All numbers are in yr-^. 



Model 


8 clocks 6 clocks 


A 
Bl 
B2 
B3 
C 


(-1.32 ± 1.46) X 10-" (-1.12 ± 1.81) x 10-" 
(-1.25 ± 1.34) X 10-" (-1.07 ± 1.60) x 10-" 
(-1.32 ± 1.40) X 10-" (-1.15 ± 1.71) X 10-" 
(-1.30 ± 1.38) X 10-" (-1.13 ± 1.67) x 10-" 
(-1.35 ± 1.46) X 10"" (-1.17 ± 1.81) x 10"" 


HBw/oD 
HBwD 

EOMS 


(-0.76 ± 0.97) X 10-" (-0.60 ± 1.06) x 10"" 
(-1.36 ± 1.44) X 10-" (-1.19 ± 1.78) x 10-" 
(-1.21 ± 1.31) X 10-" (-1.02 ± 1.54) X 10"" 


XPT+QCD 


(-1.19 ± 1.30) X 10-" (-1.00 ± 1.52) X 10-" 



Table 8: Constraints on the variation of hv once the variation of a alone is constrained from the Hg-Al 
clock. It assumes the unification relation (89). All numbers are in yr-^. 



Model 


hv 


A 
Bl 
B2 
B3 
C 


(-9.64 ±10.79) X 10-^« 
(-8.69 ±10.21) X 10-16 
(-9.27 ±10.52) X 10-16 
(-9.11 ±10.45) X 10-16 
(-9.65 ±10.71) X 10-16 


HBw/oD 
HBwD 

EOMS 


(-6.20 ±7.32) X 10-^*^ 
(-9.55 ±10.76) X 10-16 
(-8.35 ±10.01) X 10-16 


XPT+QCD 


(-8.22 ±9.91) X lO-i*" 



4.3.4 Constraint on 



hv 
Aqcd 



,a) 



As a third application, we use the fact that Aqcd = — A/^ so that the form (88) allows one to set a constraint 
on (/if /Aqcd, a) independent of any hypothesis on unification and thus does not require knowledge of 
the parameters R and S. 

Figure 3 compares the 99% C.L. constraints obtained from the combination of 6 and 8 experiments 
for each model. Again, we see that the Hg-Al experiment dominates the collective limit. 



5 Application to astrophysical systems and Discussion 

5.1 Astrophysical systems 

Several different types of observations of astrophysical systems involving quasar absorption spectra are 
subject to a similar analysis that has been applied to atomic clocks. Indeed, there are four distinct 
combinations of physical parameters which depend on g^. 

• The comparison of UV heavy element transitions with the hyperfine H I transition allows one to 
set constraints on 



X = a gpfi, 



(101) 
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Figure 3: Comparison of the 99% C.L. constraints on (/iv/Aqcd, a) for the 9 models with 8 clocks (left) 
and 6 clocks (right) 

since the optical transitions are simply proportional to Rcc- It follows that constraints on the time 
variation of x can be obtained from high resolution 21 cm spectra compared to UV lines, e.g., of 
Sill, Fell and/or Mgll. The recent detection of 21 cm and molecular hydrogen absorption lines 
in the same damped Lyman-a system at Zabs = 3.174 towards SDSS J1337+3152 constrains [87] 
the variation x to 



Ax/x = -(1.7 ±1.7) X 10" 



z 



3.174. 



(102) 



• The comparison of the Hi 21 cm hyperfine transition to the rotational transition frequencies of 
diatomic molecules allows one to set a constraint on 



y = gpa 



(103) 



The most recent constraint [88] relies on the comparison of two absorption systems determined 
both from Hi and molecular absorption. The first is a system at z = 0.6847 in the direction of 
TXS 0218+357 for which the spectra of C0(l-2), i3cO(l-2), C ^'^0(1-2), CO(2-3), HCO +(1-2) 
and HCN(l-2) are available. They concluded that 



Ay/y = (-0.16 ±0.54) x 10" 



0.6847. 



(104) 



The second system is an absorption system in the direction of PKS 1413+135 for which the molec- 
ular lines of C0(l-2), HCO +(1-2) and HCO +(2-3) have been detected. The analysis led to 



Ay/y = (-0.2 ± 0.44) x 10"^ 



0.247. 



(105) 



The ground state, ^11^/2 J = 3/2, of OH is split into two levels by A-doubling and each of these 
doubled levels is further split into two hyperfine-structure states. Thus, it has two "main" lines 
(AF = 0) and two "satellite" lines (AF = 1). Since these four lines arise from two different 
physical processes (A-doubling and hyperfine splitting), they enjoy the same Rydberg dependence 
but different gp and a dependences. By comparing the four transitions to the H I hyperfine line, 
one can set a constraint on 

F = 5p(aVA^)'-". (106) 

Using the four 18 cm OH lines from the gravitational lens at z ~ 0.765 toward PMN JO 134-0931 
and comparing the H 1 21 cm and OH absorption redshifts of the different components allowed one 
to set the constraint [89] 



AF/F = (-0.44 ± 0.36 ± 1.0; 



syst. 



X 10" 



0.765, 



(107) 
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where the second error is due to velocity offsets between OH and H I assuming a velocity disper- 
sion of 3 km/s. A similar analysis [90] in a system in the direction of PKS 1413+135 gave 



AF/F = (0.51 ± 1.26) X 10" 



0.2467. 



(108) 



The satellite OH 18 cm lines are conjugate so that the two lines have the same shape, but with 
one line in emission and the other in absorption. This behavior has recently been discovered at 
cosmological distances and it was shown [91] that a comparison between the sum and difference 
of satellite line redshifts probes the variation of 



G = 5'p(aV^) 



1.85 



(109) 



From the analysis of a system at z ~ 0.247 towards PKS 1413+135, it was concluded [92] that 
\AG/G\ = (2.2 ± 3.8) x 10"^ while a newer analysis [93] gave 



|AG/G| 



-1.18 ±0.46) X 10" 



(110) 



It was also applied to a nearby system [94], Centaurus A, to give \AG/G\ < 1.16 x 10 at 
z ~ 0.0018. 

These constraints are summarized in Table 9. 

Table 9: Constraints on the variation of different combinations of g^, ^ and a from astrophysical obser- 
vations. 



Combination 


K 


A. 


Ao 


Constraints (yr ^) 


redshift 


X = gpa'^^J. 


1 


1 


2 


-(1.7 ±1.7) x 10-^ 


3.174 


y = 5?"^ 


1 





2 


(-0.16 ±0.54) X 10-5 
(-0.2 ±0.44) X 10-'"^ 


0.6847 
0.247 


F-.9p(aVA^)i-57 


1 


-1.57 


3.14 


(-0.44 ± 0.36 ±1.0syst) X 10-5 
(0.51 ±1.26) X 10-5 


0.765 
0.2467 


G^g,{ay^^)'■^' 


1 


-1.85 


3.70 


(-1.18 ±0.46) X 10-5 
(0±1.16) X 10-5 


0.247 
0.0018 



5.2 Astrophysical constraints 

In contrast to our analysis of atomic clocks, we cannot combine the astrophysical observations because 
they have been obtained from different systems at different redshifts and at different spatial locations. 
However, as we have done previously (but without the g^ and b terms), we show in Table 10 the enhance- 
ment factor for the analysis of the 4 types of combinations of absorption spectra. We emphasize that the 
enhancement factor is always larger than unity (except for y in the EOMS and ^PT+QCD models). As 
last example of the power of coupled variations. Table 1 1 compares the constraints on the variation of 
a that can be obtained under the assumption that gp and /i are constant with the assumption of coupled 
variations based on unification. As one can see, in many cases the limits are improved by an order of 
magnitude. 



5.3 Discussion 

In this article, we have discussed the effect of a correlated variation of fundamental constants, focusing on 
the gyromagnetic factors gp and gj^. These parameters are particularly important to interpret electromag- 
netic spectra, and thus to derive constraints on the variation of fundamental constants from atomic clock 
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Table 10: Value of the parameter Ca for the 4 combinations of constants that can be constrained by 
astrophysical observations, assuming R = 30 and S = 160 (left) and value of the enhancement factor 
Cc,/Xa (right). 





X 


y 


F 


G 




X 


y 


F 


G 


A 


34.10 


15.49 


-12.60 


-17.25 


A 


17.05 


7.74 


-4.01 


-4.66 


Bl 


20.72 


2.10 


-25.98 


-30.64 


Bl 


10.36 


1.05 


-8.27 


-8.28 


B2 


24.96 


6.35 


-21.73 


-26.39 


B2 


12.48 


3.18 


-6.92 


-7.13 


B3 


23.52 


4.91 


-23.18 


-27.83 


B3 


11.76 


2.45 


-7.38 


-7.52 


C 


31.21 


12.59 


-15.49 


-20.14 


C 


15.60 


6.30 


-4.93 


-5.44 


HBw/oD 


2.46 


-15.76 


-43.23 


-47.77 


HBw/oD 


1.23 


-7.88 


-13.77 


-12.91 


HBwD 


27.00 


8.78 


-18.69 


-23.23 


HBwD 


13.50 


4.39 


-5.95 


-6.28 


EOMS 


17.63 


-0.60 


-28.07 


-32.61 


EOMS 


8.81 


-0.30 


-8.94 


-8.81 


XPT+QCD 


16.96 


-1.26 


-28.73 


-33.27 


xPT+QCD 


8.48 


-0.63 


-9.15 


-8.99 



Table 1 1 : Comparison of the constraints obtained from astrophysical systems with and without assump- 
tion on unification for model A. 



Combination 


independent (yr ^) 


correlated (yr ^) 


redshift 


X 


(-8.50 ±8.50) X 10-' 


(-4.98 ±4.98) X 10-» 


3.174 


y 


(-0.8 ± 2.7) X 10-6 


(-1.03 ±3.49) X 10-^ 


0.6847 




(-1.0 ±2.2) X 10-6 


(-1.29 ±2.84) X IQ-'^ 


0.247 


F 


(-1.40 ±3.38) X 10-6 


(3.49 ±8.44) X 10-^ 


0.765 




(1.62 ±4.01) X 10-6 


(-0.40 ±1.00) X 10-6 


0.2467 


a 


(-3.19 ±1.24) X 10-6 


(6.84 ±2.67) X 10-^ 


0.247 




(0±3.14) X 10-6 


(0±6.73)x 10-^ 


0.0018 



experiments and from quasar absorption spectra. As discussed, there is an important model-dependence 
in the computation of the gyromagnetic factors in terms of the quark masses and QCD scale. 

When applied to the interpretation of atomic clock experiments, we have shown that in general the 
constraints on the variation of a are sharper than that under the assumption that g^, g^, b and /i are 
constant, but this is not a systematic conclusion as we have exhibited models in which the variation of 
a stays the same or is even weaker due to cancellations in the sensitivity to a. The constraints on the 
variation of a should then be taken with care. In many cases, they may be stronger than reported, but 
they may be weaker as well. This points to the need to better understand the fundamental physics needed 
to calculate baryon magnetic moments. Any limit which depends on (/p „ will be subject to the type of 
uncertainties discussed here. 

Fortunately, the tightest constraint arises from the Hg-Al clock experiments, that does not depend on 
gp, gn, b or /i. As a consequence, we have been able to independently set a bound on the variation of 
hv from the combination of the other experiments. While this bound is still model-dependent, we have 
shown that it is always smaller than 



(hv) 



hv 



< 2.0 X 10" 



-15 



yr 



(111) 



for the models we have considered in this article. 

Our analysis also applies to astrophysical system and to quasar absorption spectra. We have shown 
that the enhancement factor is almost always larger than unity. 



REFERENCES 30 

Acknowledgments 

We would like to thank X. Cui, J. Ellis, M. Peskin, M. Srednicki, A. Vainshtein, and M. Voloshin for 
helpful discussions. The work of FL and KAO was supported in part by DOE grant DE-FG02-94ER- 
40823 at the University of Minnesota. JPU was partially supported by the ANR/Thales. 

References 

[1] p. Jordan, Die Naturwissenschaften 25, 513 (1937); 
M. Fierz, Helv. Phys. Acta 29, 128 (1956); 
C. Brans, and R. Dicke, Phys. Rev. 124, 925 (1961); 

R.H. Dicke, in DeWitt, CM. and DeWitt, B.S., eds.. Relativity, Groups and Topology. Relativite, 
Groupes et Topologie, Lectures delivered at Les Houches during the 1963 session of the Summer 
School of Theoretical Physics, University of Grenoble, pp. 165313, (Gordon and Breach, New 
York; London, 1964); 
J. D. Bekenstein, Phys. Rev. D 25, 1527 (1982). 

[2] J. K. Webb, V. V. Flambaum, C. W. Churchill, M. J. Drinkwater, et al, Phys. Rev. Lett. 82, 884 
(1999), [arXiv:astro-ph/9803165]; 

M. T Murphy etal, Mon. Not. Roy. Astron. Soc. 327, 1208 (2001), [arXiv:astro-ph/0012419]; 
J. K. Webb, etal, Phys. Rev. Lett. 87, 091301 (2001), [arXiv:astro-ph/0012539]; 
M. T. Murphy, J. K. Webb, V. V. Flambaum, C. W. Churchill, et al, Mon. Not. Roy. Astron. Soc. 
327, 1223 (2001), [arXiv:astro-ph/0012420]. 

[3] M. T. Murphy, J. K. Webb and V. V. Flambaum, Mon. Not. Roy. Astron. Soc. 345, 609 (2003), 
[arXiv:astro-ph/0306483]. 

[4] H. Chand, R. Srianand, P. Petitjean and B. Aracil, Astron. Astrophys. 417, 853 (2004), [arXiv:astro- 
ph/0401094]; 

R. Srianand, H. Chand, P Petitjean and B. Aracil, Phys. Rev. Lett. 92, 121302 (2004), [arXiv:astro- 
ph/0402177]. 

[5] R. Quast, D. Reimers and S. A. Levshakov, Astron. Astrophys. 415, L7 (2004), [arXiv:astro- 
ph/0311280]. 

[6] M. T Murphy, J. K. Webb and V. V. Flambaum, Phys. Rev. Lett. 99, 239001 (2007), 
[arXiv:0708.3677 [astro-ph]]. 

[7] R. Srianand, H. Chand, P Petitjean and B. Aracil, Phys. Rev. Lett. 99, 239002 (2007). 

[8] J. R. Ellis, S. Kalara, K. A. Olive and C. Wetterich, Phys. Lett. B 228, 264 (1989). 

[9] B. A. Campbell and K. A. OUve, Phys. Lett. B 345, 429 (1995), [arXiv:hep-ph/9411272]. 

[10] P Langacker, G. Segre and M. J. Strassler, Phys. Lett. B 528, 121 (2002), [arXiv:hep-ph/01 12233]; 
T. Dent and M. Fairbairn, Nucl. Phys. B 653, 256 (2003), [arXiv:hep-ph/01 12279]; X. Calmet and 
H. Fritzsch, Eur. Phys. J. C 24, 639 (2002), [arXiv:hep-ph/0112110]; 
X. Calmet and H. Fritzsch, Phys. Lett. B 540, 173 (2002), [arXiv:hep-ph/0204258]; 
T. Damour, R Piazza and G. Veneziano, Phys. Rev. Lett. 89, 081601 (2002), [arXiv:gr-qc/0204094]; 
T. Damour, F Piazza and G. Veneziano, Phys. Rev. D 66, 046007 (2002), [arXiv:hep-th/0205111]. 

[1 1] M. Dine, Y. Nir, G. Raz and T. Volansky, Phys. Rev. D 67, 015009 (2003), [arXiv:hep-ph/0209134]. 



REFERENCES 31 

[12] J. R. Ellis, K. A. Olive and Y. Santoso, New J. Phys. 4, 32 (2002), [arXiv:hep-ph/0202110]. 

[13] V.V. Dixit and M. Sher, Phys. Rev. D 37 (1988) 1097. 

[14] J.-P Uzan, Rev. Mod. Phys. 75, 403 (2003), [arXiv:hep-ph/0205340]; 
J.-P Uzan, AIP Conf. Proc. 736, 3 (2005), [astro-ph/0409424] ; 
J.-P Uzan, Space Sci. Rev. 148, 249 (2010), [arXiv:0907.3081]; 
G.F.R. EUis and J.-P Uzan, Am. J. Phys. 73, 240 (2005), [gr-qc/0305099]. 

[15] J. P Uzan, Living Rev. Rel. 14, 2 (2011), [arXiv: 1009.55 14 [astro-ph.CO]]. 

[16] K. Ichikawa and M. Kawasaki, Phys. Rev. D 65, 123511 (2002), [arXiv:hep-ph/0203006]. 

[17] C. M. Muller, G. Schafer and C. Wetterich, Phys. Rev. D 70, 083504 (2004), [arXiv:astro- 
ph/0405373]; 
T. Dent, S. Stern and C. Wetterich, Phys. Rev. D 76, 063513 (2007), [arXiv:0705.0696 [astro-ph]]. 

[18] A. Coc, N. J. Nunes, K. A. Olive, J. P Uzan, et al, Phys. Rev. D 76, 023511 (2007), [arXiv:astro- 
ph/06 10733]. 

[19] V. V. Flambaum and E. V. Shuryak, Phys. Rev. D 65, 103503 (2002), [arXiv:hep-ph/0201303]; 
V. F. Dmitriev and V. V. Flambaum, Phys. Rev. D 67, 063513 (2003), [arXiv:astro-ph/0209409]; 
V. V. Flambaum and E. V. Shuryak, Phys. Rev. D 67, 083507 (2003), [arXiv:hep-ph/0212403]; 
V. F Dmitriev, V. V. Flambaum and J. K. Webb, Phys. Rev. D 69, 063506 (2004), [arXiv:astro- 
ph/03 10892]; 

J. C. Berengut, V. V. Flambaum and V. F Dmitriev, Phys. Lett. B 683, 1 14 (2010), [arXiv:0907.2288 
[nucl-th]]. 

[20] S. J. Landau, M. E. Mosquera, C. G. Scoccola and H. Vucetich, Phys. Rev. D 78, 083527 (2008) 
[arXiv:0809.2033 [astro-ph]]. 

[21] M. K. Cheoun, T. Kajino, M. Kusakabe and G. J. Mathews, arXiv: 1104.5547 [astro-ph.CO]. 

[22] K. A. Ohve, M. Pospelov, Y. Z. Qian, A. Coc, et al, Phys. Rev. D 66, 045022 (2002), [arXiv:hep- 
ph/0205269]. 

[23] K. A. Olive, M. Pospelov, Y Z. Qian, G. Manhes, et al, Phys. Rev. D 69, 027701 (2004), 
[arXiv:astro-ph/0309252]. 

[24] T. Dent, S. Stern and C. Wetterich, Phys. Rev. D 78, 103518 (2008) [arXiv:0808.0702 [hep-ph]]; 
T. Dent, S. Stern and C. Wetterich, Phys. Rev. D 79, 083533 (2009) [arXiv:08 12.4 130 [hep-ph]]. 

[25] M. Nakashima, K. Ichikawa, R. Nagata and J. Yokoyama, JCAP 1001, 030 (2010) 
[arXiv:09 10.0742 [astro-ph.CO]]; 

C. J. A. Martins, E. Menegoni, S. GaUi, G. Mangano and A. Melchiorri, Phys. Rev. D 82, 023532 
(2010) [arXiv:1001.3418 [astro-ph.CO]]. 

[26] S. Ekstrom, A. Coc, P. Descouvemont, G. Meynet, et al, Astron. Astrophys. 514, 62 (2010), 
[arXiv:09 11.2420 [astro-ph.SR]]; 
A. Coc, et al, Mem. Soc. Astron. Ital. 80, 658 (2009). 

[27] V. V. Flambaum, [arXiv:physics/0302015]; 

V. V. Flambaum and A. F Tedesco, Phys. Rev. C 73, 055501 (2006), [arXiv:nucl-th/0601050]. 



REFERENCES 32 

[28] V.V. Flambaum, D.B. Leinweber, A.W. Thomas, and R.D. Young, Phys. Rev. D 69, 115006 (2004), 
[arXiv:hep-ph/0402098]. 

[29] S. Bize, R Laurent, M. Abgrall, H. Marion, et al, J. Phys. B: At. Mol. Opt. Phys. 38, S449 (2005), 
[http ://arXiv.org/abs/phy sics/0502 117]. 

[30] T. H. Dinh, A. Dunning, V. A. Dzuba, V. V. Flambaum, Phys. Rev. A79, 054102 (2009), 
[arXiv:0903.2090 [physics. atom-ph]]. 

[31] J. C. Berengut, V. V. Flambaum, E. M. Kava, [arXiv:1109.1893 [physics. atom-ph]]. 

[32] J.D. Pretage, R.L. Tjoelker, and L. Maleki, Phys. Rev. Lett. 74, 3511 (1995). 

[33] V.A. Dzuba, and V.V. Flambaum, Phys. Rev. A 61, 034502 (2001). 

[34] V.A. Dzuba, V.V. Flambaum, and M.V Marchenko, Phys. Rev. A 68, 022506 (2003), 
[http://arXiv.org/abs/physics/0305066]. 

[35] V.A. Dzuba, V.V. Flambaum, and J.K. Webb, Phys. Rev. A 59, 230 (1999), 
[http ://arXiv.org/abs/phy sics/980802 1 ] . 

[36] V.V. Flambaum, in Laser Spectroscopy, P. Hannaford, et aZ.Eds. (World Scientific, 2004) p. 47, 
[http://arXiv.org/abs/physics/0309107]. 

[37] M. Fischer, et al, Phys. Rev. Lett. 92, 230802 (2004), [http://arXiv.org/abs/physics/0312086]. 

[38] S. Bize, S.A. Diddams, U. Tanaka, C.E. Tanner, et al, Phys. Rev. Lett. 90, 150802 (2003), [ 
http://arxiv.org/abs/physics/02 121 09] . 

[39] T.M. Fortier, N. Ashby, J.C. Bergquist, M.J. Delaney, et al, Phys. Rev. Lett. 98, 070801 (2007). 

[40] E. Peik, B. Lipphardt, H. Schnatz, C. Tamm, et al, Proc. of the 11th Marcel Grossmann meeting, 
Berhn, 2006, [http://arXiv.org/abs/physics/0611088]. 

[41] E. Peik, B. Lipphardt, H. Schnatz, T Schneider, et al, Phys. Rev. Lett. 93, 170801 (2004), 
[http://arXiv.org/abs/physics/0402 1 32] . 

[42] S. Blatt, A.D. Ludlow, G.K. Campbell, J.W. Thomsen, et al, Phys. Rev. Lett. 100, 140801 (2008), 
[http://arXiv.org/abs/080 1 . 1 874] . 

[43] V.A. Dzuba, and V.V. Flambaum, Phys. Rev. A 81, 052515 (2010), [http://arxiv.Org/abs/1003.l 184]. 

[44] A.T Nguyen, D. Budker, S.K. Lamoreaux, and J.R. Torgerson, Phys. Rev. A 69, 022105 (2004), 
[http://arxiv.org/abs/physics/0308104]. 

[45] A. Cingoz, A. Lapierre, A.-T. Nguyen, N. Leefer, et al, Phys. Rev. Lett. 98, 040801 (2008), 
[http://arxiv.org/abs/physics/0609014]. 

[46] T. Rosenband, et al. Science 319, 1808 (2008). 

[47] R.I. Thompson, Astrophys. Lett. 16, 3 (1975). 

[48] A. Shelnikov, R.J. Butcher, C. Chardonnet, and A. Amy-Klein, Phys. Rev. Lett. 100, 150801 (2008), 
[http://arXiv.org/abs/0803.1829]. 

[49] See, for example, D. Griffiths, Introduction to elementary particles (Weinheim, Germany: Wiley- 
VCH, 2008). 



REFERENCES 33 

[50] H. Fritzsch and G. Eldahoumi, Mod. Phys. Lett. A 24, 2167 (2009), [arXiv:0906.1139 [hep-ph]]. 

[51] H. Y. Cheng, Phys. Lett. B 219, 347 (1989). 

[52] J. R. EUis, A. Fersd, and K. A. Olive, Phys. Lett. B 481, 304 (2000) [arXiv:hep-py0001005]. 

[53] J. R. Ellis, K. A. Olive, and C. Savage, Phys. Rev. D 77, 065026 (2008), [arXiv:080 1.3656 [hep- 
ph]]. 

[54] H. Hellmann, Z. Phys. 85, 180 (1933); 
R. P Feynman, Phys. Rev. 56, 340 (1939). 

[55] J. Gasser, Annals Phys. 136, 62 (1981). 

[56] R. J. Crewther, Phys. Rev. Lett. 28, 1421 (1972); 

M. S. Chanowitz and J. R. Ellis, Phys. Lett. B 40, 397 (1972); 

M. S. Chanowitz and J. R. Ellis, Phys. Rev. D 7, 2490 (1973); 

S. L. Adler, J. C. Collins, and A. Duncan, Phys. Rev. D 15, 1712 (1977); 

J. C. Collins, A. Duncan, and S. D. Joglekar, Phys. Rev. D 16, 438 (1977); 

N. K. Nielsen, Nucl. Phys. B 120 (1977) 212. 

[57] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Phys. Lett. B 78, 443 (1978). 

[58] L. F Li and T. P Cheng, [arXiv:hep-ph/9709293]. 

[59] J. Gasser and H. Leutwyler, Phys. Rept. 87, 77 (1982). 

[60] H. Leutwyler, Phys. Lett. B 378, 313 (1996), [arXiv:hep-ph/9602366]. 

[61] B. Borasoy and U. G. Meissner, Annals Phys. 254, 192 (1997), [arXiv:hep-ph/9607432]. 

[62] J. Gasser, H. Leutwyler, and M. E. Sainio, Phys. Lett. B 253, 252 (1991); 
M. Knecht, PiNNewslett. 15, 108 (1999), [arXiv:hep-ph/9912443]; 
M. E. Sainio, PiNNewslett. 16, 138 (2002), [arXiv:hep-ph/01 10413]. 

[63] K. Nakamura, et al, [Particle Data Group], J. Phys. G 37, 075021 (2010). 

[64] T. Kunihiro and T. Hatsuda, Phys. Lett. B 240, 209 (1990); 
T. Hatsuda and T. Kunihiro, Z. Phys. C 51 (1991) 49; 
T. Hatsuda and T. Kunihiro, Nucl. Phys. B 387, 715 (1992). 

[65] A. De Rujula, H. Georgi, and S. L. Glashow, Phys. Rev. D 12, 147 (1975). 

[66] A. Manohar and H. Georgi, Nucl. Phys. B 234, 189 (1984). 

[67] L. Y. Glozman and D. O. Riska, Phys. Rept. 268, 263 (1996), [arXiv:hep-ph/9505422]. 

[68] N. Isgur and G. Karl, Phys. Rev. D 18, 4187 (1978); 

N. Isgur and G. Karl, Phys. Rev. D 19, 2653 (1979) [Erratum-ibid. D 23, 817 (1981)]. 

[69] J. T. Goldman and R. W. Haymaker, Phys. Rev. D 24, 724 (1981); 
M. K. Volkov, Annals Phys. 157 (1984) 282; 
D. Ebert and H. Reinhardt, Nucl. Phys. B 271 (1986) 188. 

[70] Y Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961); 
Y Nambu and G. Jona-Lasinio, Phys. Rev. 124, 246 (1961). 



REFERENCES 34 

[71] S. Klimt, M. F. M. Lutz, U. Vogl, and W. Weise, Nucl. Phys. A 516, 429 (1990); 
U. Vogl, M. F. M. Lutz, S. Klimt, and W. Weise, Nucl. Phys. A 516, 469 (1990). 

[72] H. Georgi and A. Manohar, Phys. Lett. B 132, 183 (1983). 

[73] S. Weinberg, Physica A 96, 327 (1979); 

J. Gasser and H. Leutwyler, Annals Phys. 158, 142 (1984); 

J. Gasser and H. Leutwyler, Nucl. Phys. B 250, 465 (1985); 

J. Gasser, M. E. Sainio, and A. Svarc, Nucl. Phys. B 307, 779 (1988); 

A. Krause, Helv. Phys. Acta 63 (1990) 3. 

[74] For a pedagogical review, see, for example, S. Scherer, Adv. Nucl. Phys. 27, 277 (2003), 
[arXiv:hep-ph/0210398]. 

[75] V. Bernard, Prog. Part. Nucl. Phys. 60, 82 (2008), [arXiv:0706.0312 [hep-ph]]. 

[76] H. Georgi and L. Randall, Nucl. Phys. B 276, 241 (1986); 

H. Georgi, Phys. Lett. B 298, 187 (1993), [arXiv:hep-ph/9207278]. 

[77] E. E. Jenkins and A. V. Manohar, Phys. Lett. B 255, 558 (1991). 

[78] E. E. Jenkins, M. E. Luke, A. V. Manohar, and M. J. Savage, Phys. Lett. B 302, 482 (1993) 
[Erratum-ibid. B 388, 866 (1996)], [arXiv:hep-ph/9212226]. 

[79] U. G. Meissner and S. Steininger, Nucl. Phys. B 499, 349 (1997), [arXiv:hep-ph/9701260]. 

[80] L. Durand and P Ha, Phys. Rev. D 58, 013010 (1998), [arXiv:hep-ph/9712492]. 

[81] S. J. PugUa and M. J. Ramsey-Musolf, Phys. Rev. D 62, 034010 (2000), [arXiv:hep-ph/99 11542]. 

[82] T. Fuchs, J. Gegelia, G. Japaridze, and S. Scherer, Phys. Rev. D 68, 056005 (2003), [arXiv:hep- 
ph/0302117]. 

[83] L. S. Geng, J. Martin Camalich, L. Alvarez-Ruso, and M. J. Vicente Vacas, Phys. Rev. Lett. 101, 
222002 (2008), [arXiv:0805.1419 [hep-ph]]. 

[84] L. S. Geng, J. Martin Camalich, and M. J. Vicente Vacas, Phys. Lett. B 676, 63 (2009), 
[arXiv:0903.0779 [hep-ph]]. 

[85] D. B. Leinweber, D. H. Lu, and A. W. Thomas, Phys. Rev. D 60, 034014 (1999), [arXiv:hep- 
lat/98 10005]; 

E. J. Hackett-Jones, D. B. Leinweber and A. W. Thomas, Phys. Lett. B 489, 143 (2000), [arXiv:hep- 
lat/0004006]. 

[86] J. F. Donoghue, B. R. Holstein, and B. Borasoy, Phys. Rev. D 59, 036002 (1999) [arXiv:hep- 
ph/9804281]; 

D. B. Leinweber, A. W. Thomas, and R. D. Young, Phys. Rev. Lett. 92, 242002 (2004), [arXiv:hep- 
lat/0302020]. 

[87] R. Srianand, N. Gupta, P. Petitjean, P. Noterdaeme, et al, Mon. Not. R. Astron. Soc. 405, 1888 
(2010), [arXiv: 1002.4620]. 

[88] M.T. Murphy, J.K. Webb, VV Flambaum, M.J. Drinkwater, et al, Mon. Not. R. Astron. Soc. 327, 
1244 (2001), [astro-ph/0101519]. 



REFERENCES 35 

[89] N. Kanekar, C.L. Carilli, G.I. Langston, G. Rocha, et al, Phys. Rev. Lett. 95, 261301 (2005), 
[astro-ph/05 10760]. 

[90] J. Darling, Astrophys. J. 612, 58 (2004), [astro-ph/0405240]. 

[91] J.N. Chengalur, and N. Kanekar, Phys. Rev. Lett. 91, 241302 (2003), [astro-ph/0310764]. 

[92] N. Kanekar, J.N. Chengalur, and T. Ghosh, Phys. Rev. Lett. 93, 051302 (2004), [arXiv:astro- 
ph/0406121]. 

[93] N. Kanekar, et al, Astrophys. J. 716, L23 (2010). 

[94] N. Kanekar, Mod. Phys. Lett. A 23, 2711 (2008), [arXiv:0810.1356]. 



